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^ This paper investigates the benefits of cooperation in large wireless networks with multiple sources 

^> and relays, where the nodes form an homogeneous Poisson point process. The source nodes may 

lO dispose of their nearest neighbor from the set of inactive nodes as their relay. Although cooperation can 

potentially lead to significant improvements on the asymptotic error probability of a communication pair, 

^^ relaying causes additional interference in the network, increasing the average noise. We address the basic 

question: how should source nodes optimally balance cooperation vs. interference to guarantee reliability 

in all communication pairs. Based on the decode-and-forward (DF) scheme at the relays (which is near 

optimal when the relays are close to their corresponding sources), we derive closed-form approximations 

to the upper bounds on the error probability, averaging over all node positions. Surprisingly, in the small 
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OO node-density regime, there is an almost binary behavior that dictates -depending on network parameters- 
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ff^ the activation or not of all relay nodes. 
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I. Introduction 

In order to cope with increasing traffic demands future wireless networks must employ strate- 
gies which provide an optimal use of resources such as bandwidth and power. Over the past 
decade there has been a great interest in cooperative networks [IJ where relay nodes can 
be exploited as a means to increase throughput and reliability. Although the capacity of the 
single -relay channel [[2| remains unsolved and its optimal coding scheme unknown, significant 
progress has been made in quantifying the performance gain obtained through cooperation. 
The information-theoretic research in this topic was mainly focused on simple networks with 
few nodes or fixed topologies where perfect channel state information (CSI) is available to all 
terminal nodes. Finding explicit capacity regions of large networks may be -if feasible- very 
hard. To tackle this limitation, spatial models employing tools from stochastic geometry and 
graphs provide a comprehensive framework for the analysis of large wireless networks with 
little interference management [|J|, Q. In this setup, the interference between users is treated as 
noise [i5T| whose statistical properties depend on the particular spatial distribution of the nodes 
and fading realizations of the wireless paths. 

In this paper we study upper bounds -in terms of the average outage probability (OP)- on 
the asymptotic error probability of large-scale decentralized wireless networks in which the 
sources are aided by nearby relays. More precisely, we consider a network formed by two 
types of clusters: source-relay-destination nodes, which use a cooperative transmission scheme 
and clusters with source-destination pairs which employ simple direct transmission (DT). We 
will concentrate in the special class of networks where the relays (when available) are in the 
proximity of their corresponding sources. The rationale behind this assumption comes from 
the fact that when a source attempts to transmit, the nodes in its vicinity are the ones which, 
with a high probability, become aware of the attempts of the source and can provide help. It 
is known that when a relay is located in the vicinity of the source, the decode-and-forward 
(DF) transmission scheme [[T|, Q is capacity achieving. For this reason we will assume that 
the cooperative transmission scheme employed by a source-relay-destination triplet will be DF 
working in a full-duplex fashion. The network will be modeled as an independently marked 
homogeneous Poisson point process (HPPP), limited by the signal-to-interference ratio (SIR), 
where signal attenuation occurs both through path loss and slow fading. An outage event is 
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declared whenever the distribution of nodes and/or fading causes the target rate to be higher 
than the achievable rate. Hence the probability of these events (OP) becomes an upper bound 
on the asymptotic error probability of every pair of communicating nodes, which is the true 
performance metric of interest. 

The central motivation behind the analysis of the asymptotic error probability of large wireless 
networks in the presence of relay nodes is to provide some understanding on the limitations and 
benefits of cooperation in such decentralized scenarios. In fact, the advantage of cooperation 
for an individual source-destination link was widely studied in the past years, addressing both 
theoretical and practical issues [[T|, Q, [|7|. Although it is clear that relays can significantly 
improve the rate and reliability of a single source-destination pair, in large wireless networks 
where several source-relay-destination clusters are present the activation of relay nodes may 
drastically increase the overall interference. This occurs because as more relays are activated, 
other destinations and their corresponding relays will observe an increase in their interference 
levels. Thus, we see that although at a local level cooperation is beneficial, at a global scale it 
could be harmful. A deeper understanding of this balance between cooperation and interference 
generation is an important topic to be studied, because it could shed light on the usefulness of 
cooperation in interference-limited wireless networks of the type considered in this paper where 
the relays are in the proximity of their sources. 

Nevertheless a complete unified answer to this delicate balance is out of the scope of this 
work and is certainly very difficult to obtain, mainly because it depends heavily on the way in 
which the sources select the relays, which could be influenced by local inhomogeneities in the 
network, local network knowledge available at nodes, etc. In this paper, we aim to answer this 
question for the simple network model introduced in Section |II] with the relays positioned near 
their corresponding sources. It is assumed that CSI is not available at the transmitting nodes 
which is often the case in decentralized wireless networks without feedback. Each candidate 
relay node decides whether to be active or not in a random manner. All the nodes make their 
decisions independently of each other, and of all network parameters. That is, the decision is 
simply made through a Bernoulli experiment with success probability given by pr. 

Regarding the spatial positions of the relay nodes, two different scenarios are investigated: 

• The relays are located at a fixed relative distance with respect to their sources, 

• Every source has its nearest neighbor (NN) (among some set of inactive nodes) associated 

Febi-uary 27, 2013 DRAFT 



as a potential relay, which is obviously a reasonable assumption when using the DF scheme. 
The first case, although only of theoretical interest, is considered basically for mathematical 
tractability and because it gives some useful insights. In the second case the relays are randomly 
located, which means that in all clusters they are at different positions and hence averaging over 
all spatial positions of relays is needed to derive the OP. If a source is allowed to use its associated 
relay, cooperation will take place based on the DF scheme while if the relay is not activated, the 
source will use direct transmission (DT) to communicate with its intended destination. In this 
way the transmission scheme is a mixed one. To be fair, it should be mentioned here that there 
exist more sophisticated ways to decide whenever or not a relay node must be activated, e.g., 
by taking into account all effects mentioned above. The probabilistic model for activating relays 
that we assume is on one hand the simpler in mathematical terms to deal with, but on the other 
it still captures the balance between cooperation and interference generation via the parameter 
Pr- For instance, with < pr < 1 each source-destination pair in the network is then able to 
resign some of it local performance (choosing not to use a relay and using DT) for the greater 
good (introducing less interference into the network). In this way, by studying the effect of pr 
in the asymptotic error probability of a typical cluster, we will draw some conclusions about 
this important problem. 

A. Related Work 

Over the past years, performance gains of cooperative communications in relay networks were 
widely studied from an information-theoretic perspective. Since the seminal work of Cover-El 
Gamal [|2|, several contributions have been published on the subject. More recently, the emphasis 
was put in studying the performance of wireless relay channels where outage performance and 
ergodic rates of fading channels with Gaussian noise have been derived (see [[T|, [[6|-[[8| and the 
references therein). Among these valuable studies, the only impairments to the communication 
were due to additive Gaussian noise and fading, and very little attention was paid to the effect of 
the interference generated (or suffered) by the users. However, interference is probably the major 
impairment in wireless networks, specially in networks with little control and high mobility. 

The study of the capacity of general wireless networks taking into account the interference 
generated by the different users was pioneered by the seminal work of Gupta- Kumar [|9|, where 
the concept of transport capacity and fundamental scaling laws on the network throughput were 
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obtained considering only point-to-point coding. In [10| multiuser achievability regions were 
obtained and it was shown that for some special wireless networks significantly better scaling 
laws on the network throughput, with respect to the case in ^, are possible. Further progress 
was done in yjj where new scaling laws are derived using coherent multistage relaying with 



interference subtraction and in [12|, where extensions to fading channels were obtained. 

In all the above works, the achievability of the scaling laws derived is obtained for extremely 
regular networks which is seldom the case in practice. In that sense, stochastic geometry and 



point processes [13|, [14| are not only elegant mathematical frameworks but also seem to be 
useful tools to deal with more realistic network models, where the spatial position of nodes and 
the effect of interference can be incorporated in a probabilistic manner [|3|. Although several 
types of point processes can be used to model different kind of networks, it is the HPPP which 
has received more attention. Although other types of point processes could provide more realistic 



models [15|, the extended use of HPPP comes from the possibility to obtain simple closed form 



results in several cases of interest. The quantity called transmission capacity (TC) was introduced 



in [16 1 in order to include outage probability constraints in the scaling behavior. Several results 
have been obtained, through the use of the TC, for several practical situations, as multiple input- 
multiple output capable users in wireless networks [17], decentralized power control [[18j, etc. 
(for a review of several other important results we refer the reader to see [j4| and references 
therein). 

B. Main Contributions 

The main contributions are twofold. First, tight approximations for the OP, which lend them- 
selves to analysis, are derived. Emphasis is put on the small node-density regime (SNDR), where 
the density of active transmitters is low and which is a natural network operating point. Secondly, 
we determine the relay activation probability which provides the optimal balance between the 
reliability gain obtained from cooperation and the overall interference introduced in the network. 
From a different perspective, this addresses also the question whether or not to cooperate. The 
central conclusion we draw here is that, under quite general assumptions and for a large range 
of network parameters, the optimal activation probability pr is either or 1. This implies that 
for all clusters in the network it is best to fully cooperate or not at all. In other words, there is 
no overall gain if some clusters are permitted to cooperate and others not. 
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The paper is organized as follows. In section |II] a general and a mathematical description of 
the network model are presented. We also discuss the DF scheme and its achievable rate in the 



assumed network model. In Section III we present results corresponding to the OP performance 
when only one cluster has a relay, and the other users in the network use only DT. Although 
several results exists for the outage probability performance in fading relay channels under 
Gaussian noise [|6|, ||7|, there are no results, at least to our knowledge, where the communication 
impairments come also from a network of interferers. It will be shown that even in this simpler 
case the expression of the OP cannot be obtained in closed form and some upper bounds, valid 
when the relay is close to the source or when the outage probability is sufficiently small, are 



provided. In Section IV the more important case where there may be a relay available for 
each user in the network is treated. Several interesting situations are treated: collocated relays, 
fixed but not collocated relays, and random relays distributed with the NN distribution. Network 
parameter regions where the optimal p,. shows a binary behavior are characterized. Moreover, 
the network parameter regions where p.,. = 1 is optimal are identified. Numerical simulations are 
also presented in order to validate the derived analytical results. In section |V] a short discussion 
on the possible use of local CSI at the relays for more elaborated activation schemes is presented. 



Section VI provides some concluding remarks. Finally, several longer mathematical proofs are 



grouped together by section and deferred to the appendices. 

Notation 

We will denote M, C and M^, the real numbers, complex numbers and the real plane respec- 
tively. The canonical euclidean norm will be denoted as || ■ ||. With (•)* we denote complex 
conjugation and with 3fJ(-) the real part of complex number. Px {■} and Ex [•] denote probability 
and expectation operators with respect to the distribution of the random variable X. In a similar 
way Px|y {"1^} and Ex|y [-11^] denote probability and expectation operators with respect to the 
distribution of the random variable X conditioned to random variable Y . With A" we denote 
the complement of the set A. We shall also use the big O notation: f{x) = O {g{x)) as a; — )■ Xq 
if there exists M > and such that |/(x)| < M\g{x)\ is some neighborhood of xq. 
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II. General considerations and network model 

Consider a spatial network model in M^ in which source nodes generate messages and attempt 
to transmit them to intended destinations, either through a direct link, in which case the destina- 
tions only receive symbols from their source, or using others nodes that act as relays. Every relay 
aids a single source node and is assumed to be located in its vicinity, acting as a secondary full- 
duplex transmitter sharing the same time slots and frequency band. This setup allows the nodes 
to be grouped into clusters formed by a source-destination pair or by a source-relay-destination 
triplet, if the source has an associated relay, as shown in Fig. [T] 

The network is constructed in a two-stage fashion starting from the a set of nodes $. In the 
first stage, the source nodes are selected from the set $ and the rest is left without a specific role. 
At the end of this stage, the original set of nodes $ is divided into two disjoint sets, namely, the 
set of sources $s and the set of inactive nodes $„ which do not have a role yet. We will assume 
that the set of all nodes $ forms a homogeneous Poisson point process (HPPP) of intensity A 



1 14 1, and that the role in the first stage is assigned through independent marking. This means 
that at the end of the first stage, the set of nodes $ can be thought as the superposition of 
two independent HPPPs, namely ($3, $»„), with intensities (A^, Aj„), such that A = A^ + Aj„. It 
should be clear that this can also be thought as that the medium access control (MAC) protocol 
working in the modeled network is a slotted ALOHA [3], such that at each slot every node 
decides whether to transmit or not. It is assumed that the intended destinations for the sources 
nodes are not part of the process <l>. 

At the second stage every source node will be associated or not to an active relay node in 
its vicinity. Hence the problems to be dealt with are how such decision should be made and 
how a relay should be chosen among a set of suitable candidates. It is clear that an adequate 
decision rule for relay activation should take into account local criteria, to satisfy the clusters 
needs, and global ones, to achieve certain goals such as network throughput. As a matter of 
fact, if all relay nodes had full CSI of all (including destination) nodes involved in the network, 
they could construct a "global" consensus to evaluate (depending on fading, node proximities, 
etc.) the best way to associate relays to source-destination pairs and at the same time to avoid 
causing too much interference to each other. Nevertheless, this does not apply to our setting since 
we shall assume only little -local- information is available at relays, no information about the 
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Source 



Fig. 1. General additive relay channel, gsr, Prd, and Qsd are the channel gains between source and relay, relay and destination 
and source and destination, respectively. The interferers i, j and k and their corresponding channel gains toward the relay and 
the destination are also shown. 



destination nodes and no CSI is available at the source nodes. Although these assumptions may 
appear to be restrictive, they are quite realistic in decentralized wireless networks. For instance, 
we assume that the relays solely decide their activation randomly and with the same probability 
and independently of all network parameter^ From a mathematical standpoint, this means that 
each node of the set of sources $s is assigned with a Bernoulli mark of success probability Pr, 
indicating whether the relay associated to a source is activated or noj^ It will be assumed that 
the corresponding sources remain oblivious of the relay decisions. Regarding the relay selection, 
we shall consider two situations: 

• For a source node at position x, the relay is positioned at a fixed relative position from the 
source denoted hy r = x + k, where A; is a fixed two-dimensional vector. 

• The relay is chosen by the source from the set of unassigned nodes $j„. 

The first case is considered for mathematical tractability and because some useful insights 
about the general problem can be obtained. With respect to the second case, it is reasonable to 
assume that relays will choose their sources using the scarce local knowledge of the network 

'in section V we briefly consider the case in which the relays use local CSI in order to decide their activation. 
^Note that this model can also handle the possibility that a source does not have a potential relay available. 
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they possess. The relay assignment is a very complex issue, since on one hand a one to one 
correspondence between the sources and their relays has to be established, and on the other hand 
this correspondence must be established so that the pairs are appropriately located for effective 
cooperation. Since the focus will be on networks where the relays are in a vicinity of their 
sources, we propose that each source will attempt to associate its nearest neighbor (NN) from 
the set $«„. This assignment rule does not avoid the possibility of two sources being associated 
to the same node as a relay since the same unas signed node could be the closest neighbor of one 
or more sources. Nevertheless, assuming a reasonably small node-density functioning regime the 
density of active sources (As) will be much smaller than that of the potential relays (Xm) and 
the probability of conflicts will be small. 

With these simplifications, the existence and location of all the relays become independent for 
each source and we can simply include the position of the NN of each source as an independent 
mark distributed with the NN probability distribution. In this way, we do not have to deal with 
conflicts and a possible dependence between the relay assignments in different source nodes, 
simplifying the model considerably. 

As a final remark, it is important to note that in each cluster we have a mixed scenario, i.e. 
with probability pr a cooperative scheme is used while DT will be implemented with probability 
(1 — Pr)- Notice that the relay activation probability pr becomes the crucial parameter which 
dictates the balance between cooperation and the aggregate interference in the network. As Pr 
increases a greater proportion of clusters enjoy the benefit of cooperation, supporting a greater 
rate or having a better quality of service. But at the same time more interference is introduced in 
the network. Therefore, the main question would be: which is the optimal activation probability 
at a global scale? From a practical point of view, this question is relevant because its answer 
would increase our understanding of the limits of cooperation in interference-limited networks. 
A priori, the answer is not obvious and in the rest of this paper we will try to give some insights 
about it, at least for our simplified model. 

A. Mathematical Description 

Based on our simplifying hypothesis described above, we now provide a more rigorous 
description of the network model. The network is modeled as an independently marked HPPP: 

^ s ^ \_\^i) ^Xii Ki, hxiri iT-Xidi i^kif, iT-kid) f ■ (t) 
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satisfying the following requirements: 

• The positions of the sources constitute the HPPP $s = {xi} of intensity A^. 

• The Bernoulli random variable (RV) e^.. indicates that source Xi uses a relay, according to 
an event with probability pr. 

• The relay associated with the source at Xj is located at Xj + ki, where fcj is a degenerate 
RV or with the NN from $,„ distribution. All source-destination distances are D. 

• All nodes transmit with unit power while the power received at y by a transmitter at x is 
\hxy\'^l'xy where l^y = \\x — y||~" (a > 2) is the usual path loss function and |/ixyP is the 
power gain of Rayleigh fading] with unit mean. This is equivalent to saying that h^y are 



complex, circular [19|, zero-mean Gaussian RVs. 
• An additional source with the same marks as the others, independent of the point process 
$s and with its destination at d = [D,0], is added at the origin. The position of the 
relay for this source node will be r (distributed as a degenerate RV or distributed with the 
NN distribution). The coefficients j/i^rP^ l^rdP and |/isdP model the source-relay, relay- 
destination, and source-destination fading coefficients of this cluster, respectively. This 



cluster is independent of the point process $5. Slyvniak's Theorem [13|, [14| guarantees 
that the study of this cluster's behavior will be representative of the behavior of any other 
similar cluster in the network and hence it can be considered as a "typical cluster". 
• hx^r, hkiv and hxid, hk^d model the fading gains between each source and its relay and the 
relay of the source at the origin, and between each source and its relay and the destination 
of the source at the origin, respectively. Thus the gains shown in Fig. [T] correspond to 

|„ 12 I/, 127 |„ 12 I/, 12/ „„J |„ 12 \l, 127 

\ysr\ \i''sr\ ''sri |i/sd| \i''sd\ ''sd "iiU \yrd\ \'''rd\ ''rd- 

It will be assumed that the HPPP $s remains fixed during the transmission time. 

Remark 2.1: The path loss Ixy function is a very accurate model of wireless power propagation 
when the distance between the points x and y is sufficiently large (far- field approximation). 
However, it is not a good model when x is close to y because of the existence of a singularity 



which makes it unbounded. A more suitable choice would be pO| |: 

1 + ||,T — y\\°' 

^In this paper we focus on Rayleigh fading but, the results can be extended to the case of Nakagami fading |17|. 
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which avoids the singularity and has the desired behavior when ||x — y|| is large. In despite of 
this, the above function would increase the analytical complexity with small benefits when the 
density of active transmitters As is small [|4), [21 1 (typical operation regime). 

We shall assume that at each position y in the space a background complex, circular and 
zero-mean Gaussian noise Zy is present. This noise has variance 1/SNR and it is statistically 
uncorrelated in time. There is no need in specifying any spatial correlation on such background 
noise. Then, by conditioning on the marked HPPP, the signals received at the relay and destination 
-associated with the source node at the origin- can be written as: 



Yr 



f^sr^s 



\r\\ 2 



+ 



s 



iT'Xir^x 



•^i 



^ h \\Xi + ki 



r 2 



(3) 



Y. 






r 2 



r 



d\[- 



i-XiG^s 



h Y 

T* ' Ci 'T' ■ 

\xi - d\\'' 



+ ^x,: 



hkidXk, 



JU'i ~P fX/ii Qj 



+Za 



(4) 






where, for shortness, we have dropped the dependence of the signals on the messages to be 
transmitted and the discrete time indices for the block codewords. We have denoted with (X^, Xj) 
the complex, circular and zero-mean Gaussian signals at the source and relay (if any is associated) 
for the cluster at the origin, and (Xj,^, X^J the corresponding signals for the other clusters in the 
network. Source and relay nodes in each cluster are assumed to not perform any joint decoding 
of their own messages with the messages of other clusters in the network. Therefore, they will 
simple take those interference signals as pure noise. Therefore, the relay and the destination 
associated with the source at the origin see aggregate interferences Zr and Z^, respectively, 
yielding the following lemma. 

Lemma 2.1: Let a > 2, then for almost all realizations of the point process $s, the aggregate 
interferences Z^- and Z^i are zero-mean complex circular Gaussian variables whose conditional 
variances are respectively given by 



E 

iiXiG^s 

E 



I'^Xirl 



\Xj — r\ 



^Xid\ 



\xi - d\\' 



+^x, 



+ea 



\h 



kir\ 



2^{K,rhl^p} 



Ixi + fci-rll" 



\h |2 



\xi — r\\ 2 \\xi + ki — r\\ 2 



\xi + ki- dW"" \\xi - d\\ 2 \\xi + ki- d\\ 2 



(5) 
(6) 
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Proof: The proof follows from the fact that when a > 2, I^ and /^ are finite for almost 



every realization of $s . This can be shown using the Laplace functional of $s 1 13 1, [22 1, and 
the following functions 

x/ , r. r. ^ ^ ^ l^ll' f I^Sp 23? {/li/l*p} \ ^^^ 

/ X, 5,/;;, /ii, /i2, hs, K) = — + e j — — r- +- -^ — — -^ (7) 

U; _ ^ a yiix-^-k — rW^ \\x — r\\2\\x + k — r\ 



2 



and g{x,e,k,hi,h2,h3,h4) defined in an similar form. As the signals (X^.,Xfc.) and (X^^ ,X, 






are Gaussian and independent when i j^ j, from equation (13) below it can be shown that the 
partial sums (through a proper enumeration of the points of the particular realization of $s) in 
Zr and Zd are Gaussian with variances given by the corresponding partial sums in Id and Ir- 
Thanks to the finiteness of Id and J^, this family of Gaussian distributions have the tightness 
property (Theorem 25.10 in [23]) which allow us to obtain the desired result. ■ 

B. Problem Statement and Bounds on the Asymptotic Error Probability 

Our goal will be to study the asymptotic error probability performance of cooperative com- 
munication schemes, which could make use of the presence of potential relay nodes. In order to 
do that let us now define the notion of code which is the same for any cluster in the network, 
similarly to that of point-to-point codes in [|5|: 

Definition 2.1 (single-relay code): A single-relay code-C„(n, R, Mn) for a set {1, ... , M„} of 
uniformly distributed messages W consists of a set of randomly and independently generated 
complex Gaussian codewords X", each according to n i.i.d. draws of a Gaussian RV of unit 
variance, a decoder mapping W : C" i — > {1, . . . , M„} U {^}, and a sequence of relay mappings 
ft ■ C*^^ I — > C constrained to produce i.i.d. complex Gaussian RVs of unit variance, for 
t = {1,. . . ,n}. The smallest asymptotic average (over all random parameters) probability of 
error of, for example, the source-destination pair at the origin and which will be representative 
of any source-destination pair in the network, is given by 

Pe{R) = inf I limsupP^"^(Vr y^ W\Cn) 

where 6 condense all the randomness in the model: 



lim inf — log M„ > R 

n— 5>oo n 



Q = \^s,hsr,hrd,hsd,r,eo\ . (8) 
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Notice that the source is unaware of the instantaneous interference, path loss attenuation and 
fading coefficients involved, which implies that the error probability cannot be made arbitrary 
small with the code-length. 

The asymptotic error probability respect to a source node located at the origin can be upper 
bounded using any code C„ as follows: 



PJR) < inf 



Pe {0{R)} + limsupP^"^ \ W ^ W\Cn, 0%R) 



(9) 



where cr(0) is the cr-algebra generated by and 0{R) denotes any outage event. By choosing 
adequately the outage event 0{R) for a given code C„, it turns out that the second term on the 
right-hand side of (|9]) can be made arbitrary small. That is, for any e > 0: 

limsupP^j"^ Iw ^ W\Cn,0^{R)\ < e. (10) 

In this way, fixed a rate R, the asymptotic error probability Pe{R) is dominated by the outage 
probability (OP) Pe {0{R)} of the corresponding achievable rate. The OP is a useful perfor- 
mance metric which was extensively employed to characterize performance in a Poisson field 
of interferers, jointly with the associated metric of transmission capacity (see Q, [ 16 1, [22| and 
the references therein). 

The average error probability of any single -relay code-Cn(ri, i?, M„) defined as above can be 
lower bounded as follows [|24l: 



P,{R) > miFe{R>mm{I{Xs;YrYd\Xr),I{XsXr;Yd)}} 

>miFe{R>IiXsXr;Yd)}. (11) 

The rationale for the second inequality comes from the fact that for the class of networks 
considered here, where the relays are in a vicinity of the sources, the limiting term in the cut-set 
bound is /(X^X^; Yd). 

In our context of Poisson networks we shall also consider the scaling behavior of the error 
probability with the density of interferers. We have the following definition: 

Definition 2.2 (small node-density regime): We define the small node-density regime (SND) 
as the situation in which XsirD"^ — )• characterized by the following metric: 

/.(i?)^ liin ^^. (12) 
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As XgirD^ is the average number of active nodes in circle of radius D, k(R) gives the asymptotic 
behavior of the smallest probability of error for a code with rate R, when the number of active 
users in disk of radius D centered at the destination goes to zero. In other words, k(R) measures 
the reliability of the communication when the closest point to the destination is the source that 
intends to transmit to it, and the destination does not see the interference from other nearer nodes 
which would have stronger channels because of their reduced path loss. In this way k(R) is a 



metric, similar to the standard notion of diversity gain [25| used for fading Gaussian channels. 



that measure the robustness against the interference generated by the network. 

C. Achievable Bounds on the Asymptotic Error Probability 

The main coding strategies for relay networks were introduced in the seminal work by Cover- 
El Gamal Q. There have been three dominant relaying paradigms: decode-and- forward (DF), 
compress-and-forward (CF), and amplify-and-forward (AF). Decode- and-Forward (DF) scheme 
allows the relay to decode the messages sent by the source, re-encode them, and forward them to 
the destination. In the special case where the memoryless relay channel is physically degraded 
the achievable rate using DF is in fact the capacity. However, in the general case, it turns out 
that, depending on the channel parameters, there is no a unique scheme maximizing the rate. In 
general, roughly speaking, if the source-relay channel has less noise than the source-destination 
one, the relay can perfectly decode the messages intended to the destination without introducing 
any bottleneck in the information flow. This does not hold true for all cases since the imposition 
of full decoding at the relay could be a strong one. Although some variants as partial-decode- 
and- forward (P-DF) [2| may ameliorate this inconvenience, we will focus on the DF scheme. 
The potential improvement of P-DF essentially relies on a careful optimization of the code at 
the encoder, which cannot be correctly done in our setting due to the lack of CSI at the source. 
Moreover, our interest will be in networks where the relays are located in some vicinity of 
their corresponding sources, generally closer to them than to the associated destinations. In that 
situation, as path loss has a major role in channel quality, assuming that the distance source-to- 
relay is smaller in average than that of source-to-destination, the DF scheme should not introduce 
significant penalties with respect to other schemes. 

There exist several coding and decoding strategies based on block-Markov superposition 
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encoding that achieve the same DF rate. In [2], cooperative strategies use irregular encoding, 
random binning and successive decoding at the destination. However, this method is not well- 
suitable for the present scenario. Helpfully, the same rate can be derived with regular encoding 
and sliding-window decoding p6| at the destination. Another cooperative scheme that permits 



to achieve same rate is regular encoding and backward decoding p7| , p8| . We assume that 
transmissions are done based on regular encoding and sliding-window decoding [[T| with Gaussian 
signaling. The ri-length random codewords at each source and its associated relay is written as 

X:{w,_,,w,) = ^{1 - \p\^)Xl\wi) + pX;\w,.,), X,"(t^,_i) = X2"(u;,„i), (13) 

for messages Wj G {1, . . . , 2"^^^} with wo = wb+i = 1 and each block z = {1, . . . , B}, where 
Xi and X2 are independent complex, circular Gaussian RVs with unit variance and p is the 
correlation coefficient between source and relay signals Xg and Xj.. Then, conditioned on a 
particular realization of $s, we can guarantee that for the relay channel associated with the 
source located at the origin the following rate is achievable using the DF scheme [|2|: 

D • Lf\^sr\'^hr{l-\P\'^)\ ^(\hsd?lsd+\Kd\\d + '^\/hJ7d^{phsdKd)\\ 

Rdf = max mm < C -, , C -, > , 

P6c,|p|<i \ V /r + l/SNR )' V /d + l/SNR jy 

(14) 
where C{u) = log2(l + u). 
Define the outage event when the relay is present as Oi{R) = {A{R, p) U B{R, p)} with 

A{R,p) = {\hsr\^lsr (1 - IpI') < T(/, + 1/SNR)} , 

BiR,p) = \\hjnsd+\hrd\Xd + 2Vu7d'iiiphsdK^)<T{h + l/SNR)Y 



for which condition pO] ) holds true. The event A{R, p) means that the relay is in outage while 
B{R, p) means that the destination is in outage while source and relay cooperate. T = 2^ — 1 
and p is the complex correlation coefficient between the transmitted source and relay symbols. 
The error probability when the source at the origin has an associated relay is bounded by: 

Pe(i?ko = 1) < Pout,DF(i?) = inf Pe{Oi(i?)ko = l} 

pec,|p|<i 

= inf PeM(i?,p)Ui3(i?,p)|£o = l}. (15) 

pec,|p|<i 

We also define the outage event Oq{R) = AoriR) for the case in which there is no relay and 
thus the source simply uses DT with 
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for which again condition (fTOl) holds true. For the source at the origin and without relay, the 



error probability is upper bounded by 

Pe{R\e^ = 0) < Pout,DT(i?) = Fe {C»o(i?)ko = 0} = Pe {ADT{R)\eo = 0} , (17) 

and for the general case with mixed situations, the upper bound writes as 

Pe{R) < Pout,m:x(i?) = inf inf p,Pe {Oi{R)\eo = l} + (l-p.)P© {Oo{R)\€o = 0} . (18) 

o<Pr<ipec,|p|<i 

Remark 2.2: It is important to mention that the events A{R, p) and B{R, p) are the outage 
events only if regular encoding with sliding-window decoding or regular encoding with backward 
decoding are used. If irregular encoding with successive decoding is used we should consider 
another event regarding the estimation by the destination of the binning index sent by the relay, 
taking into account the rate of that transmission which is different from R. 

Remark 2.3: The DF scheme with backward or sliding-window decoding at the destination 
are oblivious [8] to the potential absence of the relay. That is, in the absence of the relay (and 
without the source knowledge of that event), the DF scheme does not degrade the performance 
with respect to DT In this way, the coding scheme employed by the source can always be the 
same (equation ([T3])) independently of the activation or not of the associated relay. 

Remark 2.4: The rate Rdf does not depend on the correlation between the noises or interfer- 
ences at relay and destination. This is true because the correlation between received signals at 
the relay and the destination becomes irrelevant when full decoding at the relay is imposed. As 
a matter of fact, this is not the case for the CF and AF schemes where the correlation between 



the noises can increase or decrease the corresponding achievable rate |29| 



III. A SINGLE RELAY CHANNEL AND MULTIPLE SOURCES 

In this section we will consider the case when only the source at the origin has a relay 
associated with probability one and the other sources in the network transmit without the help 
of a relay. The goal of this section is to provide a quantitative analysis of how cooperation 
improves performance in a single cluster and how interference reduces its benefits. In the case in 
which only Gaussian background noise and Rayleigh fading is considered (without interference) 
very interesting gains have been observed in terms of OP [[6| [[7|. The setup with interference 
in which only one source has a relay can be considered as a best case scenario from the point 
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of view of this source, and it will give insight to the general problem of DF performance in 
interference-limited networks. Although this case may seem like a direct extrapolation of the 
case with background noise, in this scenario we must average over all possible configurations 
of interfering nodes, considering numerous situations in which communications are severely 
impaired due to the presence of heavy interference. As we will see, this results in performance 
gains which are not as large as in the case in which only fading and noise are considered. 

In the following we shall consider that there is no background Gaussian noise (SNR —^ oo), 
and that the only impairment is the interference generated by the other nodes in the network. 
The motivation for this is twofold. On one hand, it permits us to focus on the more important 
effect of the aggregate interference, which could be several times greater than the background 



noise [22|. On the other hand, not considering the background noise will lead to more compact 
expressions which are easier to analyze. Nevertheless, the results obtained could be modified 
without too much effort in order to include the effect of background Gaussian noise. 

We start by deriving the OP in the case in which the relay of the cluster at the origin is at 



a fixed known location (Theorem (3.1 )). Although a closed form expression can be derived in 
terms of the joint Laplace transforms of the interferences at the relay and at the destination, 
it involves certain integrals which cannot be evaluated in closed form [30]. Therefore we find 



upper bounds on the OP which are tight when the OP is small (Theorem 3.2| ). In addition, we 



study the optimal correlation coefficient p between the symbols transmitted by the source and the 
relay, and find that in the SND regime its value is zero. Afterwards we study the case in which 
this single relay is randomly distributed according to the nearest neighbor distribution of the 
source. Direct evaluation of the OP in this case can not be obtained through exact expressions, 
but starting from the bounds previously obtained for the fixed relay we develop upper bounds 
for the OP in this scenario, which are tight when the OP is small (Theorem 3.3| ). This allows 



for a characterization of the OP and a comparison to DT in a more realistic scenario. 
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A. Relay at a Fixed Position 

The outage probability when the source at the origin is the only one having a fixed relay at 
a known position r can be written aqj 

Pout,DF(i?) = inf {1 - Pe {A'iR, p) n B\R, p)}} . (19) 

P6C,|p|<l 

In the following, in order to simplify the expressions, we will drop the inf operation and it will 
be assumed implicitly. From the definitions of A{R, p) and B{R^ p) we can therefore write: 

TO tr>\ ^ la \ \U \2 -^ '^^r \hsd\ hd + IKdl kd + '^Vhdlrd^iphgdKd) ^^1 

[ isr (1 - |pr) h J 

(20) 
It is interesting to mention that the OP has two different expressions according to the relay 
position r and the correlation coefficient p of the symbols transmitted by the source and the 
relay. The following Theorem, which deals with the expression of the OP in this setup, includes 
both expressions. However, as we shall see in what follows, working with only one of them is 
enough for characterizing the OP behavior. 

Theorem 3.1 (OP for a single fixed relay): The outage probability for the DF protocol when 
only the cluster at the origin has a relay and | |r — (i| | 7^ Z^ or p 7^ is: 

Pou,DF(i?) = 1 ^^£/„/,(T/p2, T/ps) + —^^^Cj,j^{T/pu T/ps), (21) 

P2 — Pi P2 — Pi 



where: 



and: 



{hd + Ird) — {{hd — Ud) + ^lsdlrd\p\ ) ,^^. 

/^l = ^ , (22) 

1 /9 

_ {hd + Ud) + {{hd — Ird) + ^lsdlrd\p\ ) ,r,r,. 

/^2 — , K^J) 

P3 = /.. (1 - IpP) , (24) 



Ci^j^ {u,,u;2) := E^^ [e-i'^^^^+'^^^r)^ , u,, U2 G C, (25) 



''since in this case only the source at the origin has a relay associated at position r with probability one, the randomness in 
the model is given by O = \^b, hsr, hrd, hsd \- Similarly, the marks in "l>s associated with the relays are not considered, and 
the interference at the destination associated with the source at the origin can be written as: 

"^ ^ lis, -dlh' 

with a similar expression for the interference at the relay's position. 
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with 'Si.{ui} ,'Si.{uJi} > is the joint Laplace transform of the interference at the relay and at 
the destination. In addition, when ||r — (i|| = D and p = 0, we have pi = ^2 and the outage 
probability is: 

T dCj^j^{ui,T/ij,3] 



IPout,DF(-R) 

Proof: Defining: 



1 + 



/ii 



dui 



-^h,ir{T/f^i,T/fi3). 



uJi=T/fj.i 



V '■= \hsd\ hd + \hrd\ Ird + '^y hdKd^iphsdKd) 1 

and considering that hgr and V are independent of each other and of $s we have that: 

,2 ^ TIr V 



P, 



out,DF 



{R) 



E, 



1-E, 



P 



hsry\^s 



P 



hsr\^s 



3 -L 
TL 



\hj^ > 



f^3 T 

<I,Ip 



/^3 



v\^, 



Ml 


7 


{^ 


>Id 


q 



1-E 



*. 



e ''^^Fv{Th) 



(26) 
(27) 

(28) 

(29) 
(30) 



where -FV(") is the complementary cumulative distribution function (CCDF) of V and yU3 is given 



by (24). Notice that in general V has the distribution of the sum of two independent exponential 
RVs with different means. The exception to this is when p = and ||(i — r|| = D. In that 
case, the means of the exponential RVs are the same so V follows a Gamma distribution with 



2 degrees of freedom. In precise terms, in Appendix D-A we show that the CCDF is: 



At2e 



-"/M2- 



/^le 



-"/mi 



Fv{u) 



A«2-A«i 

(l + M//^i)e-"/^i pi 



Pi 7^ P2 

/^2 



(31) 



where pi and p2 come from ([22]) and (23). Straightforward calculations show that pi = p2 only 



if p = and ||(i — r|| = D. Now, replacing that expression in (30) and using the definition of 



the Laplace transform (25) we obtain (21). To obtain (26) we replace (31) in (30) and use the 
fact that: 



(iE^^ \e~ 



(wi/tj+W2/r)l 



duji 



E^^ [-/, 



dC 



-(aJl7'd+aJ2-fr) 



]■ 



(32) 



Remark 3.1: Notice that the OP in this scenario is a split function because the fading RV 
V given by ( p7| ) has the same property. However, the expression of the OP when p = and 
||(i — r|| = D, given by (26), can be obtained by continuously extending the other expression 



February 27, 2013 



DRAFT 



20 

pT] ) at these points. Therefore, in what follows we shall focus our interest on ( pT[ ) which fully 

characterizes the OP. 

Now we consider the case of the simplified path loss function: 

Corollary 3.1: For the simplified path loss function, when \\d — r\\ ^ D or p ^ the OP is: 



Pout,DF(i?) = 1 ^^e" 

/i2 -/il 



_/^2 ^-A,[5(/.-2/°+/.3-2/")+/{T/M2,r/M3)] 

where 



/i2 - yUl 



/(a;i,a;2) = / 7 n ttt^t n vr^dx, (34) 

c^^r(^)r(i-i). (35) 

S = CT^/", (36) 
and r(2;) = J^ t^^^e^^dt is the usual Gamma function. 



Proof: Use Lemma |A.1| in Appendix |A] to evaluate the Laplace transform in the previous 
Theorem. ■ 

Remark 3.2: The OP depends only on the absolute value of p and not on its phase. This is a 
consequence of the uniform phase of the Rayleigh fading coefficients. 



Notice that /(wi, U2) in (34) accounts for the statistical dependence between the interferences at 
two different locations. If in fact this two interferences were independent, then the joint Laplace 
transform would result in the product of their individual transforms and this cross-term would 
not appear. Unfortunately this term does not have a closed form and is very difficult to bound 
it tightly in a general setup for all (101,102). Therefore, although we can find numerically the 
value of f (001,002) without difficulty, we cannot provide a closed form expression for direct 
computation. Since in general we assume that the relay is in the vicinity of the source, the union 
bound on the outage probability will be tight. That is: 

Pout,DF(i?)< inf (^e{A(R,p)} + f>e{B(R,p)]), (37) 

pec,|p|<i 

will be a good approximation, since a close inspection of A(R, p) and B(R, p) shows that in 
this setting the event B(R, p) will be dominant and that A(R, p) will have a relatively small 
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probability of occurrence. The following Theorem deals with evaluating the union bound ( p7| ) 
and with the small node density regime: 

Theorem 3.2 (Upper bounds on the OP for a single fixed relay and the SND regime): When 
only the source at the origin has a fixed relay and \\r — d\\ ^ D or p ^ 0, considering the 
simplified path loss function, the OP can be upper bounded as: 



Pout,DF(i?) < inf 

pec,|p|<i 



fJ'2 - jJ'l 



(38) 



with /ii and /i2 given by ( [22| ) and ( |23] ). In addition in the small node density regime we can 
bound k{R) as: 



5 I \\r\ 

k(R) < inf - — — 



1 /^2 



1-- 



1-- 



P 



2\2/a £,2 



/i2 -/il 



Proof: For the first part it is easy to show that: 



Pe {A{R, p)} = l-Ci^ (r/(l - IpH) = 1 - e (- 



As-illrl 



2W^ 



IpP) 



(39) 



(40) 



which represents the OP of a DT from the source to the relay p| with the exception of the 
correlation coefficient p which does not introduce any difficulty. In addition, the proof that: 

1^2^-1^ (T/fi2) - /ii£/, (T//ii) 



Fe{B{R,p)} = l 



1 - 



P2 - Pi 

- r -2/a , r -2/o 

yU2 -yUl 



(41) 



(42) 



follows along the same lines as Theorem 3.1 The Laplace transform has the same expression 
as ( [40l ) (it can be found by setting (jJ2 = in ( |116| ) from lemma A.l in Appendix [A]). 

The proof of the second part follows directly from the fact that e^" = \ — u + 0{u^) and the 



definition of k{R) (12) 



Notice again, that the case where ||r — c?|| = D and p = can be treated via continuity arguments 
as mentioned above. 

Now, it is important to adjust the value of p which determines the correlation between the 
source and relay signals. We do so in the SND regime. 

Lemma 3.1 (Optimal p in the SND regime for the union bound): In the small node density 



regime, the right side of (37) is minimized by taking p = 0. 



Proof: See appendix B-A 
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Fig. 2. OP when only the source at the origin has a relay Fig. 3. Outage probability as a function of \p\ for various 
located at r. Exact expressions for DF come from j33| l while relay position and for \s — 10^'^ and Aa — 10~*. d — (10 0), 
upper bounds come from (|43|. d — (10, 0), R — 0.5, a — A. a = 4 and R = 0.5. 



Although the optimality of p = is proved for the SND regime, we will see through numerical 
simulations that it is in fact true for practical values of the OP pT) . 

Remark 3.3: As pointed out in [jT|| (see remark 42) and [|7j, using p = simplifies the 
implementation of DF because symbol synchronization between the source and its corresponding 
relay, which could be difficult to obtain in ad hoc wireless network, is not strictly required. 
In what remains of this section we focus on the case p = 0. 

Corollary 3.2: Using the simplified path loss function and when ||(i — r|| ^ D the optimal 
OP (p = 0) can be upper bounded as: 

£)Og-As(5||r-d|p _ IIt- - 



Pout,DF(i?) < 1 - e 



,-^sS\\r\\' 



1 - 



\a^-XsSD^ 



D' 



In the small node-density regime we have: 



k{R) < 



6D 



TT 



D2 



+ \\r-d\\ 



,D 



a-2 



\r — d 



|a-2 



D^ 



dt 



(43) 



(44) 



Now we present a few figures to show the behavior of the expressions derived and to compare 
the performance of DF in this case. In Fig. |2] we can see the comparison of DF versus DT, 



both through the exact numerical evaluation of the OP using (33) and with the upper bounds 



given by (43) for different relay positions, taking d = (10 0), a = 4 and R = 0.5. We can see 
that these bounds are very accurate when the OP is small as proposed. In addition, for a fixed 



February 27, 2013 



DRAFT 



23 



source-relay distance the OP increases as the relay grows further away from the destination. 
This is because the probability of B{R, p) increases as this happens. In Fig. [3] we can observe 



how the variation of the true outage probability of equation (33) as a function of \p\ for various 
values of r, using the same parameters as in figure [2} Two sets of curves are presented. One for 
the case of A^ = 10^"^, in which the OP is small, and the other for A^ = 10^^, in which the OP 
is larger. In both cases we see that p = is the optimal choice. 

B. Relay at a Random Position 

In this section we consider that only the node at the origin has a relay with probability one, 
and that this relay is chosen as the NN of the source at the origin. Thus, the expressions for the 
OP found in the previous section can be interpreted as being conditioned on the position of the 
relay so now for a random relay located at r the OP can be expressed a^ 

Pout,DF(i?)= inf Er [Fe\r {Aip, R) U Bip, R)\r}] . (45) 

pec,|p|<i 

In the previous section the OP for a fixed known location of the relay was calculated and it was 
discussed that for the small node-density regime the optimal correlation coefficient between the 
symbols of the source and the relay is in fact p = 0. Assuming the same operating condition, the 
probability inside the expectation should be minimized for all values of r when p = 0. Therefore 
in this scenario we can write: 

Pout,DF(i?) = Er [Fe\r {^(i?) U B{R) \r}] , (46) 

where we define A{R) = A{R, 0) and B{R) = B{R, 0). It is clear that the OP for the simplified 



path loss function with ||r — (i|| j^ D (and p = 0) can be calculated by averaging (33) from 



corollary 3.1 with respect to the relay position r: 



KutMR) = 1 - E, 



D' 



D" - ||r-(i||' 



-Xs[5i\\r\\^ + \\r-dn-finr-dr,T\\r\n] 



-Xs[5{\\r\\HD^)-fiTD",T\\r\n 



(47) 



D"' -\\r-d\ 
The values of r for which ||r — (i|| = D do not affect the value of the expectation since at 



those points (47) has an avoidable singularity. From the last section we know that the function 



^In this case the randomness in the model is given by 9 = •^ ^s,hsr,hrd,hgd,r> where r is now random, and we also 
remove the marks in "I>s related to the relays of the other clusters. 
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/(■, ■) cannot be computed in closed form. Moreover, it is very difficult to upper and lower 
bound accurately in a general scenario. This precludes in turn a closed form computation of the 
expectation with respect to r as it is required in this scenario. Notice also that each term inside 



the expectation does not have finite expectation separately so the expectation in (47) should be 
evaluated jointly. For this reason, to derive closed form expressions we have to find a suitable 
upper bound which can be averaged with respect to r. In the previous section we showed that 
in the SND regime the union bound is tight, and therefore it will also be useful in this scenario. 
Therefore we upper bound the OP as: 

Pout,DF(i?) < K [Peir {ARM + Peir mR)\r}] . (48) 

The next lemma deals with the nearest neighbor distribution of a HPPP which we require for 
averaging: 

Lemma 3.2: The position of the node closest to the origin of an homogeneous Poisson 
point process of intensity Au„ is a two-dimensional Gaussian r.v. with zero mean, independent 
components and variance (2Aj„7r)~\ that is, its density is: 

Mr) = ;^e"^ = A,„e-"--ll^ll'. (49) 



m 



Proof: The proof follows from the fact that the distance from the origin to the closest point 
in the Poisson process is distributed as a Rayleigh RV. Also, from the homogeneity of the point 
process, it follows that the angle of vector r should be a uniform RV in [0, 2tt), independent of 
its norm. ■ 



The expression inside the expectation in ( |48] ) is in fact ( [43] ), which we now want to average 
with respect to r. The first term is due to A{R)\r and its average has a closed form expression. 
The second term is due to B{R)\r and cannot be averaged in closed form. Moreover, it has an 
avoidable singularity at \\r — d\\ = D which further complicates matters. Therefore we propose 
to find a simpler upper bound on the probability of this event. We restrict ourselves to the case 
a >4. 

Lemma 3.3: When a > 4, Pe {B{R)\r} can be upper bounded as: 

Pe {l3iR)\r} < 1 - {N - M\\r - d\\) e'^'^^', (50) 
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where: 



A^ = 1 + K^D"^ + 



M 



\s6D [ 1 



2{K6D^f 
a 

1 - K6D^) 



a 



(51) 
(52) 



for any set of values {Xs,T,D). 



Proof: See appendix |B-B[ ■ 

The usefulness of the previous lemma relies on the fact that the bound is tight when | |r — (i| | ~ D 
(for example near the origin) and that it is linear in ||r — c?|| so it can be averaged in closed 
form. Notice that this bound becomes negative if | |r — (i| | is too large; however it is easy to show 
that if XgSD^ < 1 (even well out of the SND regime) this bound will be positive at least when 
I |r — (i| I < 2D. Since the relay will be close to the source, the contribution of the negative parts 
will be negligible after averaging. With this bound we can upper bound the OP in the small 
outage regime: 

Theorem 3.3 (Upper bounds on the OP for a single random relay): When « > 4 the OP for 
the DF protocol when only the source at the origin has a relay following the NN distribution 
and for the simplified path loss function, can be upper bounded as: 



Pout,DF(i?)<l-(iV-ME,[||r 



+ 



vrAirt + As(5' 



where A^ and M come from (51 ) and (52). In addition in the SND regime, we have: 



k{R) < - 

For the expectation we can use: 

Er[||r 
Erillr 



2<T'tn^(^_2\^r\\\r 



D^ 



a 



D 



= crinQ2,o{D/crin,0), 
< D 



'CTin 
7 1^^) (^r 



D 



(53) 

(54) 

(55) 
(56) 



where Q2,o is the (2,0) Nuttall Q-function f3T|. In particular: 



Q2,ois,0) 



\-'^{{s' + 2)lo{j]+s'h 



where Iq and Ji are the modified Bessel functions of the first kind of orders and 1. Similarly 
7(s) is defined by: 



7(s) 



1 + 



TT 



2 erf 



\V2s)_ 



(57) 
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where erf(-) is the standard error function. 



Proof: The first two terms in (53) come from using lemma 3.3 to bound Peir {'B(-R)|r}, 



which is the second term of ( |43) ), and averaging with respect to r. The third term in ( |53) ) 
follows from averaging Pei^ {A{R)\r}, which is the first term of (43), with respect to the relay 



distribution (49 ). The SND regime result once more follows from the fact that e " = l—u+0{u 



and the definition of k{R) ( 12 ). Finally for the proof of the results on the expectation of | |r — (i| | 



see appendix B-C 



Remark 3.4: ( |55| ) will provide tighter expressions than ( [56] ). However, ( [56] ) will allow further 
analysis in the next section when we shall include the possibility of a helping relay to all source 
nodes in the network. 7(5) is monotonically increasing so its minimum value is attained as 
s — 7- when 7(5) — t- (^ — l) ~ 0.3425. In addition, if s < 1/3 it is almost constant. Notice that 



if an upper bound is determined for s then 7(5) can be upper bounded by a constant so (56) 
can simplified to a linear function of a^. 

To finish this section we present a few figures to study the performance of the bounds defined 
in Theorem (3.3 ). In figure |4] we can see the OP for a random relay for different values of a^ 



obtained both through numerical simulation of (47) and by using (55) in (53) against DT In all 



cases we take a = A, R = 0.5 and d = (10,0). We observe that the bounds derived predict the 
OP in the small outage regime very accurately. It is clear that as cri„ grows, that is, as the NN 
distribution becomes less concentrated around the source, the probability of finding the relay 
close to the source decreases and on average the performance of DF will decrease. In figure 
[5] we plot the maximum achievable rates attainable through DF relative to the rate of DT for 



an OP constraint, as a function of the NN dispersion. The DF rates are obtained by using (55) 



in (53); since this is an upper bound on the OP, the plotted gains are actually a (tight) lower 
bound of the actual gains. It is clear that as the spread of the NN decreases {am/ D decreases) 
the average gains become larger. 



IV. Multiple sources and relays 

In the previous section we considered the case in which only the cluster at the origin has 
the possibility of cooperating. This was done to isolate and analyze how a situation of strong 
interference could affect the performance of the DF scheme of a given user in the network. 
Although this analysis is important in itself it only provides a local understanding of cooperation 
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Fig. 4. OP when only the source at the origin has random Fig. 5. Maximum DF rate relative to the same rate of DT 
relay. The exact OP for DF comes from numerically evaluating for a fixed OP. DF rates are obtained by using l|55|l in (|53l. 



1 47), while bounds come from using ([55| in (|53|. d = (10, 0), A = 10"". d = (10, 0), a = 4. 
R = 0.5, a = 4. 



since only one user is allowed to use a relay. This is similar to the standard AWGN fading relay 
channel. In this section we shall consider the more interesting situation in which relays are 
available for all source nodes. As we shall see, many of the tools developed in the previous 
section will aid us in what follows. We start by discussing a useful approximation which will 
lead to a mathematically tractable problem. 



A. A Useful Approximation 

Even with the simple model presented in section |ll| the problem at hand is mathematically 
difficult to attack. The expression of the OP for the transmission scheme detailed in section [ll] 



is given by equation (18). In more precise terms, that equation can be expressed as: 

inf Pn„,n,ix(i?,Pr), 



P, 



out,mix 



{R) 



0<Pr<l 



(58) 



where Pout,mix(-R,Pr) denotes the outage probability of the scheme for a fixed value of p^- Using 
the tools presented in appendix [A| it can be shown that Pout.mixl-R^Pr) can be written as: 



P, 



out,mix 



{R,Pr 



inf Pj-E, 
pec,|p|<i 



i^U,ir {T/fi2, T/ns) + \. Ci^,ir {T/f^i, T/n3 



/i2-/Xl 



/^2-/ii 

+ {l~Pr)CjAT/lsd), (59) 
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where /ii, /i2 and /is are defined in (22), (23) and ( [24] ) respectively. The interferences /^ and /^ 



are given by (Bl) and (|6]) and the two-dimensional Laplace transform Ci^j^{uJi,U2) are given by 

Cj^j^{Ui,UJ2) = exp {-XsPrt{uJi,U2,r,d, p)} exp ^-X,{1 - Pr) C{ujI^" + ul^") + f{uJi,UJ2) |, 



(60) 



where f {001,002) is given by (34) and t{ui,uj2,r, p) is given by: 

t(ui,U2,r,d,p) = / Ek[l - s{ui,x,k,d,p) s{u2,x,k,r,p)]dx. (61) 

The expectation is with respect to the distributioii|j of k and where s{oo, x, k, d, p) is given by: 

s(co,x,k,d,p) = n 77, ir, 77, r. r-7^^ — ^n 77, — ir, 77, — • (62) 

For s{uj,x,k,r, p) a similar expression holds interchanging d with r. Notice the complexity 
of the true OP. It is clear that expression ( [59] ) is not amenable for mathematical analysis. 
Moreover, obtaining ([59]) numerically implies significant computational resources because of 
the need of performing integration routines in high dimensional spaces (for example M^, if 
the relay are distributed with the NN distribution). For this reason, we shall make another 
simplifying assumption. The variances of the interferences (conditioned on $s) in the relay 
and the destination associated with the source at the origin are given by ([5]) and ([6]). These 
expressions prove to be very difficult to work with, specially if we are interested in simple final 
closed form results, which could lead to mathematically tractable analyzes. Besides the statistical 
dependence between I^ and Ir, the problem is given by the different path loss attenuations, that 
for example, the destination or the relay of the source at the origin, see from a given active node 
and its associated relay. This difference in path loss attenuation lead to integrals on M^ or M^ 
which cannot solved in closed form or be tightly bounded (from below and above). As we are 
interested in simple expressions which could shed light on the problem at hand, we shall avoid 
this problem by considering a simplified path loss attenuation: 

^ \\xi + Tki - r\\° 

^According to our model in section II k is distributed as a degenerate RV in R^ or witli tlie NN from <I>i„ distribution. 
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with r G [0, 1]. The main difference between these expressions and the ones in ^ and ^ is 
given by the fact that we have approximated the two path losses appearing in the interferences 
by a single one. For the case of the destination of the source at the origin this becomes: 

11 11 
^ ^ (65) 

Iki-rfll" \\xi + Tki-d\\''' \\xi + h - dW'' \\xi + rh - dW"" 

Similar expressions apply for the relay. Consider the term ||xj + rki — d]]^"', for what it is easy 
to show that: 

1 1 oL „ "^^",„ V (66) 



\\Xi + Th - d\\" ||x, -rf||" V \\xi-d\\ 

So, in the case \\xi — d\\ ^ Ci\\ki\\ the approximation should be good. Of course, this cannot be 
guaranteed in an exact way, because of the random locations of the nodes Xj. However, if A^ is 
small enough, this can be assured with high probability. It is easy to show that: 



P0{||s-d|| <a||A;|||} = 1 - E^ 



-AsTrallfcl 



(67) 



where the latter expectation is with respect to a generic mark k. Notice that under the assumptions 
in section |n| the distribution of k (as a fixed vector or as the NN distribution of $j„) is 
independent of the particular x at which is associated, that is, the marked HPPP assumed for 



the model is stationary [13|. When the relay position is fixed and if As||A;p is small enough, the 
probability of having a network configuration in which the approximation was not good should 
be small. In the case when k is distributed with the nearest neighbor distribution of $j„, 

Pe{||x-rf||<a||A;|||} = -^^, (68) 

and the same would be true depending on A^ and Ai„. Typically, in order to have useful levels of 
OP performance, the number of active nodes per unit area, A^, should smaller than the number 
of inactive nodes per unit of area Aj„. In that sense PeJUx — (i|| <a;||A;|||} should be small. We 
should mention, that the above reasoning is meant to be an intuitive argument and not a formal 



proof to validate (63) and (64). However, the numerical results will show that this is indeed a 
good approximation. We will use that approximation below in order to analyze the OP of the 
typical cluster and its dependence with pr. First we will analyze the problem when the relay 
is collocated with the source and then when it is not collocated but it is at a fixed position. 
Finally we shall treat the problem when the relay is distributed around each source with the NN 
distribution of $j„. In the following we shall also impose that p = and although we will be 
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working with ([63]) and ([64|) we will still use the symbols /^ and Id just to simplify the notation. 



Remark 4.1: We should note that the choice of p = will not improve necessarily the global 
performance of the network. This is because if p 7^ the Gaussian interference generated by a 
source-relay pair will be correlated. For another node subject to this interference that correlation 
could be beneficial. On the other hand, from a local point of view, p = could improve the OP 



performance of the DF scheme (see section III I so there is another interesting balance between 



local and global network performance. In this work, however, we will not consider this problem. 

B. Relays Collocated With their Sources 

We start by analyzing the case in which the relay is collocated with the source [[8| in all 
the clusters, because this setup is mathematically tractable and it will give us some interesting 
insights on the more general problem. That is, we assume that the source-relay channels in all 
the clusters are not subject to amplitude fading and only to a fixed common gain given by (7 G C 
(with the possible presence of phase fading) known to the relays. The relays are not impaired by 
the interference /^ and only by complex, circular Gaussian noise with unit variance. This could 
model the situation in which the relays have out-of-band reception and in-band transmission fy2\. 



1 33 1, that is, the source-relay links are orthogonal to the source-destination and relay-destination 
hnks, which could happen for example, if the source-relay links are wired. The relays are also 
at relative distance r (ki = r Vi) from their corresponding source nodes, in such a way that the 



approximation (64) at the destinations is valid (obviously approximation (63) is not needed in 
this case), and their channels to their corresponding destinations are subject to fading and path 
loss. As usual the sources and relays radios transmit with unit power. The sources should also 
invest some power in the transmission to the relays through their corresponding links, which is 
already considered in the channel gain g. 

It is clear that with this simplified model we still capture the effect we are interested in. 
Although the relays reception is not impaired by the network interference, in order to obtain 
local cooperative gains through DF, the relays have to transmit (in a wireless fashion) introducing 
interference in the network and degrading the global performance. 

Let us consider the Bernoulli RV Eq indicating that the source at the origin uses a relay. To 
find the OP we can condition on this RV; when the node at the origin does not use a relay then 
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the OP is the same as DT, and when it does, it is the OP of the DF protocol. So we can write 
the outage probability (with p = 0) a^ 

Fout,mix(i?,Pr) = Pei.o {AiR) U B(R)\eo = l}P.o(^o = 1)+P0|eo {ADTiR)\eo = 0}P.o(^o = 0), 

(69) 
where we used A{R) = A{R,0) and B{R) = B{R,0). The event A{R) has probability one or 
zero depending on log (1 + Ig]"^) > R or log (1 + Ig]"^) < R. As we are interested in exploiting 
the relay if it is present, we restrict ourselves to the case log (1 + \g\'^) > R. In this way, the 



relay, if present, is always able to decode the source message. Equation (69) becomes: 



Pout,mix(i?,Pr-) = Pei.o {ADT{R)\eo = 0}F,,{eo = 0) + Pe|.o {B{R)\eo = IjP.ol^o = !)• (70) 

We have the following Theorem: 

Theorem 4.1 (OP for collocated relays in all clusters): The OP when the relays are collo- 



cated with their sources in all clusters (eq. (70)) is: 

Pout,mix(i?,Pr) = Pr (1 - -^/^ {T/hd)) + (1 " Pr) 

i)" - ||r-(i||" ' D" - ||r-(i|| 



1 7~^i-^ II _7 1 1 i-ii T~^r\/ II JlWrM I ' ^ ' 



where: 



CjM) = exp j-A.Cu;^/" (^ + ^) } ' ^^^^ 



and C is defined as in p5\ and 6 is given by (36). 



Proof: The proof of the expressions in terms of the Laplace transforms can be obtained 



along the same lines as Theorem 3.1 The expression of the Laplace transform can be found in 
lemma A. 2 in appendix |A[ ■ 



Notice that the Pr appearing in the exponentials of (71 ) accounts for the additional interference 
introduced by the other active relays in the network. On the other hand, the values of pr and 1 —pr 
that weigh these exponentials account for the possible local improvements due to cooperation 
in the typical cluster. It is easy to verify that when pr = (there are no active relays with 
probability one), Pout,imx(i?, 0) = Pout,DT(^)- 



^We are making an abuse of notation in saying that l|69| is tiie true outage probability of l |59| >. This not the case basically 
for the path loss approximation mentioned in subsection |IV- A[ However, as l |69[ l will be an excellent approximation to the true 
outage probability and in order to avoid introducing a new symbol we will use the same notation for this approximation. 
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The following lemma deals with the optimal relay activation probability: 

Theorem 4.2 (Concavity of pr for collocated relays): If the following inequality holds true: 



then Pout(-R,Pr) (eq. (70)) is concave function in pr- 



Proof: See appendix C-A 



We have then the following corollary: 

Corollary 4.1 (Optimal p,. for collocated relays): For the collocated relay case introduced above, 
if the destination and relay positions satisfy: 

\s5max{\\r-df,D^} <0.38, (74) 

we have that Fout{R,Pr) attains its minimum value when pr = or pr = I. 

Proof: It follows from the fact that given a concave function /i(x) in a bounded and closed 
interval [xi,X2], its minimum is attained at xi or X2 (see Theorem 32.1 in [34]). ■ 



This means that under condition (74) the best OP performance for any cluster in the network 
can be attained when all or none of the sources decide to use their associated relays. In one 
case all the clusters will be using DF and in the other case all of them will be using DT This is 
a somewhat surprising result in the sense that in terms of the OP the best performance can be 
obtained either by full cooperation or by not cooperating at all. There is no "optimal" density of 
used relays in the network or optimal mixed behavior in the sense that some clusters would enjoy 
the advantages of cooperation while others use DT in order to balance the generated interference. 
The conditions under which pr = I or pr = would be the preferred choice depend on several 
network parameters (a, d, r, Xg, R, etc). To see that, in the following subsection we treat the 
special case when r = 0, that is when the relays are physically collocated with the sources. 



Remark 4.2: The result of lemma 4.2 is only a sufficient condition for Pout,mix(-R, Pr) to be 



concave. It is not a necessary condition. In fact, our numerical results suggest the condition 

Xs6miix{\\r-df,D^}<^, (75) 



is also sufficient for the concavity of Pout(-R,Pr)- Notice that (75) covers almost all the important 
conditions of operation. Consider that when XgSD^ = 0.5, the outage probability for DT is 
Pout,DT(-R) ~ 0.4 which already corresponds to a very poor network performance. 
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C. Relays Physically Collocated With their Sources 

In order to apply the results derived above we will consider the ideal situation in which the 
relays are physically collocated with their sources. It is easy to see that this corresponds to the 
case of a MISO system (2 transmitters- 1 receiver), and it achieves the second term in the DF 



rate expression (14 1. The same situation happened in the case treated in the last subsection. The 
difference is that both antennas did not present the same average gain due to the different path 
losses they were subject to. In this case, since both of them are located in the same position, 
they have the same average gain. However, in the cases treated before and in this one, there 
is a fundamental difference with respect to the usual MISO case: in some sources the second 
antenna might not be used. In fact, it will be used with probability pr. Therefore, the balance 
between the improved local performance and interference generation is considered through the 
use of the second transmitting antenna. The OP for this case can be exactly computed along the 



same lines of Theorem 4.1 and for that reason we skip the proof of the following result: 

Theorem 4.3 (OP for physically collocated relays): The outage probability Pout,mix(-R, Pr) for 
the case in which the relays are physically collocated with the sources is: 



■''^out.mix I, ^ 1 Pr ) ^ 






exp <! -A,(5D2 ( 1 + ^ ) \. (76) 



The following lemma identifies the exact optimal values of Pr as a function of the other 
parameters of the network: 

Theorem 4.4 (Optimal Pr for physically collocated relays): The optimal value p°P^ for the case 
that the relay is physically collocated with its source is: 

1 0<v<jf^ 



^opt 



2 



ri 2/1 



«^- Tfw<V<2 (77) 



Pr 

[ r/> 2 

where r] = XgSD"^. In the SND regime, where the optimal p^ is one: 



<R) < - f 1 - 4, ) . (78) 



TT \ a"^ 



Proof: The proof follows standard procedures of function optimization and for that reason 
is omitted. ■ 
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Fig. 6. PoLit.raix(i?,Pr) and Pout.DT(-R) as function of p,,. D ■ 
10, A, = 10"*, a = 4, R = 0.5. 



Fig. 7. Pout.mix(^,Pr) and Pout.DT(-R) as function of p^- D 
70, A^ = 10"*, Q = 4, i? = 0.5. 



Notice that in this ideal case the dependence of the optimal pr on the network parameters is 
condensed on a and r]. As Pout,DT(-R) = 1 — e"'' all the above results can be put as functions 
that depend solely on the path loss exponent and the performance of DT with the actual network 
parameters. 

In Fig. |6} |7] and [8} we see the outage probability of the mixed scheme and DT versus pr. It is 
clear that DT does not depend on pr. In Fig. [6] we see the case where the network is operating 
in < 1] < Y^ and for that reason pj. = 1 is optimal. In Fig. 17 the operating point is such 
that YzW < ^ < 2. In this case the optimal pr is not or 1 (we see that Pout mix (-R,Pr) is not 
concave). Finally, in Fig. 3 we have i] > 2 and DT is optimal network-wide. In this case, we 
also have that Pout,mix(-R5Pr) is not concave. So, in this particular case of physically collocated 
relays, in the concavity region obtained in Theorem |4.2[ pr = I is optimal. In other cases this 



is not true as we shall see later. It is important to mention that the cases in Fig. |7} and [8] occurs 
when the outage probability performance is far from practical and typical values. That is, from 
a practical standpoint, p^ = I is optimal. 

In Fig. |9] we see how the optimal pr behaves as a function of Pout.DT(-R) for various values of 
a. We can observe that, for a large range the optimal probability is equal to 1. We also see that, 
for all values of a, the region where 1 or is not the optimal value of Pr is sufficiently small. 
Finally, independently of a, when Pout,DT(-R) > 0.865 the optimal value is zero. We see there is 
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Fig. 9. Optimal relay activation probability pr as a function 
of Pout,DT(i?). D = 10, A, = 10"*, i? = 0.5. 



an almost binary behavior, but this occurs far from the typical small node-density regime. This 
also shows what was mentioned above: from a practical standpoint, for this ideal case, Pr = 1 
is optimal and cooperation is beneficial globally. 



D. Relays at a Fixed Position but not Collocated 

We shall consider now the more important case where the relay is at a fixed relative position 
with respect to its source for all the network clusters but it is not collocated. That is, the source- 
relay link, when it is used, is subject to fading, path loss and network interference. The outage 
probability is given by (|69]) (rewritten below for easy reference): 



P, 



out, mix \ 



{R,Pr) = Pei.o {A{R) U B{R)\eo = l}P.o(^o = 1)+P0|eo {AdARMso = 0}P,„(£o = 0) 



(79) 
Again, we shall use the union bound in order to provide some tractable expression. The exact 
expression of Pout,mix(-R, Pr) is given by (59). The use of that expression will complicate the 



analysis without providing any significant improvement in the accuracy. As mentioned in section 
Tin ^f '^he relay is positioned in the vicinity of its source, the union bound should be a reasonable 
and accurate expression, because the event B{R) will be dominant with respect to A{R). 
Remember also that the assumption of the relays located in the vicinity of its sources was 
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also considered in order to use the approximation in the path loss attenuation in the subsection 



IV-A In this way we can write: 



Pout,mix(i?,P.) < (Peieo {AiR)\eo = 1} + Pei.o {B{R)\eo = 1}) P(£o = 1) 

+ Pe|eoMDT(i?)ko = 0}P(£o = 0). (80) 

We have the following Theorem: 

Theorem 4.5 (OP for fixed non-collocated relays): For the situation in which the relays are 
fixed at relative position r with respect to the sources, and subject to fading and network 
interference we have: 



Pn 



c(i?,Pr) < {l-Pr) [l-Cl, (T/lsd)] +Pr 



D-CjAT/lrd) ^\\r 



'^i. (T/U) 



D^-llr-rfll" D^-llr-rfll" 

+ Pr[l-CjAT/lsr)], (81) 



where: 



£/^(wi) = exp <^ -A^Cw/" 1 + 



2/q 



2pr 



a 



(82) 



and 5 is given by (36) and C is defined as in (35). 



Proof: The proof follows the lines of Theorem 4.1 The new term Peieo {-^(-R)ko = 1}^ 



which represents the outage event at the relay, can be found following the same steps used to 
obtain Pei^o {ADTiR)\so = 1}. ■ 

It is important to notice that the term which considers the outage at the relay: 

2pr 



Pr 1 -exp<^ -Xs5\\r 



a 



(83) 



is not a concave function of pr- In this way, we cannot be sure that the upper bound to 



Pout,mix(-R,Pr) in (81) is concave in pr. However, our numerical results indicate that at least 
when As(5max {||r — rfp, D'^} < 0.38, that seems to be true. In fact, it is easy to see that for all 

Pr- 



0<Pr{ l-exp<^-A,(5||rf 1 + 



2pr 



a 



< 1 -exp<^ -As(5||r|p 1 + 



2pr 



a 



(84) 



that is, the upper bound in (81) to Pout,mix(-R,Pr) can be tightly bounded between 2 concave 
functions of Pr- Notice that the tightness is given by the quality of the link source-relay. As 



long as As5||r|p is small enough, the bounds in (84) will be very tight and the upper bound. 
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and it is very likely that Pout.mix(-R,Pr) is a concave function. Again, notice that As5||rp small 



is coherent with the use of the path loss approximation in IV-A and the union bound. In this 



way, the results of corollary 4.1 should also apply to this case under the conditions mentioned 
above. 

In order to obtain more insightful results let us introduce: 

A{pr) = ^ (l + ^) • (85) 

In general we shall be interested in A(l) = 5 (l + -), so unless we specify the value of pr we 
shall simply write A = A(l). Now lets assume that A^A = A^^ (l + ^) is small. As mentioned 
before, this will be in general the typical regime of operation. Usually ||r|| < D and \\r — d\\ 
will be comparable to D. These observations will allow us to have the following result: 
Lemma 4.1 (SND regime for fixed but non-collocated relays): In the SND regime we have: 



k{R) < m{pr), (86) 



where: 



m(pr) = — M^ - D'^-^ ^ p 



2() /WrP D"^ — \\r — dW' > 

^" Ml' II nQ-2 ^ II' "II \ J2 



r 



Proof: Follows from the fact that e"" = 1 — u + 0{v?). ■ 

In the small node density regime we can use the function m(pr) to study regions in which pr = 1 



is optimal. For that reason we have the following lemma proved in appendix |C-B[ 

Lemma 4.2 (Region of optimality of pr = 1 for fixed relays): The function m{pr) presents the 
following properties: 

• Vpr G [0,1] it is concave if the relay relative position r is contained in: 

A{a,d) = [jLeR^:\\rf<D^^—^, \\r - df = xD^\ . (88) 

• The value pr = I is optimal if the relay relative position r is contained in: 

1-x 2 



B(a,d) = [j Ir eR^:\\rf<D^' 



X2 a + 2 



d\\' = xD'}. (89) 



It is interesting to observe that B{a, d) C A(a, d), which implies that in a practical setting, 
the set of r G M^ where pr- = 1 is optimal is contained in the region where m(pr) is concave. 
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Fig. 10. Outage probability Pout.mix (^,Pr) as a function of Fig. 11. Outage probability Pout,mix(^,Pr) as a function of 
Pr when r = (1, 1) and r = (3, 3). The other parameters pr when r = (5, 4) and r — (—3, —2). The other parameters 
are d = (10, 0), A^ = 10"^, R = 0.5. are d = (10, 0), A^ = 10~"*, R = 0.5. 



We also see that lima-s>oo B{a, d) = A{a, d) implying that when the path loss exponent is large 



enough the two sets coincide. 



Remark 4.3: It can be shown that the unions in the ([88]) and ( [89| ) can be restricted to 
corresponding compact sets in M^ and not all x > 0. It can also be shown that both regions 
B{a,d), A{a,d) are compact subset in M?. 



In Fig. [To] and 11 we see the outage probability Pout,mix(-R, Pr) versus pr for different settings 
and using the numerically calculated exact outage probability (]59]) and the approximate one 



in (81). In all cases it is important to mention how tight is the approximate expression. We 
see that in some cases, the approximate expression in ( JST] ) is not a strict upper bound to 
the true Pout,mix(-R, Pr) in (59). This is because of the path loss approximation discussed in 



subsection ]IV-A In any case the curves are tight enough to be useful approximations to the true 
Pout.mixl-R^Pr)- From the figures we can also see that when the relay is sufficiently far from the 
source, the performance of the mixed scheme degrades with respect to DT. The same happens 
when the relay is located near the source but "behind" it. In the first case, the reason from 
the degradation is from the fact that the relays are far from the sources and cannot decode 
the message reliably. In this way, interference is introduced in the network due to the relays 
presence, but without enjoying any benefit from cooperation. In the second case, the relays are 
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Fig. 12. Regions of concavity and optimality of p,. 
using in{pr) in l|87| for a = 2.5 and a — 6. 



1 Fig. 13. Regions of optimality of p,. = 1 calculated using 
the exact outage probability l |59[ ), the approximate expression 

10"*, a = 4, 



I 81 1 and function m{jpr). d — (10, 0), As 
7? = 0.5. 



able to decode reliably, but the cooperative gain in the second term of (14) is not sufficient to 
deal with the aggregate interference generated by the relays. This is because, the quality of the 
link relay-destination is not good enough to provide a boost in performance due to cooperation. 



In Fig. 12 we see the regions where m{pr) is concave and where Pr = 1 is optimal for that 
function, for two different values of a. We see that that as a is increased, the source is displaced 
from the center of the regions and is positioned in a corner. In the same manner, the regions 
are expanded towards the destination. This figure clearly shows that the best relay positions 
(where pr = 1 is optimal) are biased in the direction toward the destination. Positioning the 
relays "behind" the source is not the best, specially when the path loss exponent a is large 



enough. This result agrees with the intuition about the problem. In Fig. 13 we see the regions of 



optimality of pr = I obtained using the numerical exact outage probability (59), the approximate 



expression (81) and function m{pr). We see the excellent agreement of the obtained simplified 
expressions with the true behavior given by the exact expressions. 
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E. Relays Randomly Located Around Their Sources 

In this section we assume that the potential relays are distributed according to the NN 
distribution around each source. This means that in each cluster the relay will be located at 
a different position. To derive the OP we have to average with respect to sq as we did in 
previous sections, and also with respect to the position of the NN of the origin, r. As we have 
done so far, we shall consider that p = 0, and we shall use the union bound to obtain expressions 
amenable for analysis. Using these hypothesis the OP can be written as: 

Pout,mix(i?) < P.o {^0 = 1} (E. [Pe|,o,r {A{R)\eQ = 1}] + E, [Pe|.o,. {B{R)\eo = 1}]) 

+ PeMDT(i?)ko = 0}P,o{£o = 0}, (90) 

where we have used that r is independent of sq, and that Adt is independent of r. The following 
Theorem deals with the expression of the OP and the SND regime: 



Theorem 4.6 (OP for random relays and the SND regime): The OP when the relays are se- 
lected as the nearest neighbors of each source can be upper bounded as: 

AsA(pr) 



,(i?)<(l-p. 



1 



-\s/^{Pr)D^ 



+ Pr{l 



TlXin + \s^{Pr, 



l + A,A(p,)D" 1 + 



9 / 2 — a^ 

aU 



-0 {{XsA{pr)Dy) 



as X,A{pr)D'^ = \,5 (l + ^)D^ ^ 0. In the SND regime we have: 

2\ EJ 






1 - P. 1 



JJ2 



a 



D 



(91) 



(92) 



The expectation of ||r — (i|| is given by (55) or can be upper bounded by (56). 



Proof: For the first part see appendix |C-C and the second part as usual follows from the 
fact that e" = l + u + 0(u2). ■ 

An important issue in this setup is to study the optimal relay activation probability, since now 
each relay is randomly located. In the previous section the optimal choice of pr depended on 
the set of parameters (a, d, r. As, R). In the random relay setup, the position of the relay is 
different in each cluster. Hence, the optimal relay activation probability will depend on the set 
of parameters (a, d, ain, Xs, R), that is, it will now depend on the average distance to the 
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NN with respect to the transmission distance D. In what follows we shall neglect the term 
O {{\sA(pr)D'^y) in (91) which is a natural approach when the OP is small. 



Theorem 4.7 (Concavity of the OP for random relays): When the OP is small, neglecting the 
term O {{\s^{pr)D'^Y) in (91), for each value of XgSD"^ there exists an interval < crj„ < 



ain,c{D, Xs6D'^, a) such that the OP is a concave function of pr G [0, 1]. 

In addition if \s5D'^ < 0.38 then a lower bound of aj„ ^(-0, XgSD"^, a) is obtained by solving 
the equation: 

e^^^^' {2as^ + (« - 2)7(s)s) - {XsAD^ + (4 + X,6D^)X.AD^ - (2 + 3X^60^) = 0, (93) 

where s = ain/D. Finally a closed form lower bound of crj„ c is given by: 

(a -2) 



cr,;, 



> 



(A,AD2)2 _ (4 + XjD^)XsAD^ + (2 + 3Xs6D^] 



2a 



-e 2 



XsAD^ 



8a 



(94) 



Proof: See appendix C-D ■ 

Finally the following Theorem deals with sufficient conditions for p^ = 1 to be the optimal 
relay activation probability. 

Theorem 4.8 (Optimality region of pr = 1 for random relays): When the OP is small, neglect- 



ing the term O {{XsA{pr)D'^y) in (91), for each value of XgSD"^ there exists an interval < 
o'in < o"m,i(-D, Xs5D'^, a) such that the optimal relay activation probability is pr = 1. 

In addition if XgSD^ < 0.38 then a lower bound of ain,i{D, Xs5D^, a) is obtained by solving 
the equation: 

2a {a + 2) e^^^^'s' + {a^ - 4) ^{s)s + 4 ('^XJD^ - l) = 0, 

where s = ain/D. Finally a closed form lower bound of o-jn.max is given by: 



(95) 



/ 2(l-iVD2) .,^^^. _ (a-2) ,^^. 



a{a + 2) 



10a 



(96) 



Proof: See appendix C-E 



To finish this section we present a few figures to study the expressions we just derived. In 
Fig. [14] the OP with respect to pr is plotted for two different values of a in, one in which Pr = 1 



is optimal and another one for which p^ = is the optimal point. The approximations of the 



OP come from using (55) and (56) in (91 ) and they are compared with Montecarlo simulations 
obtained by averaging 8 x 10^ realizations of the PPP using the true interferences (|5| and ([6]), 
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taking d = (10, 0), A^ = 10^^, R = 0.5. We see that the approximations derived with the 
simplified interferences are in excellent agreement with the actual OP derived with the more 
complex interferences. 



In Fig. 15 we study the maximum values of (Jm/ D for the concavity of the OP and for Pr = 1 



to be optimal in the SND regime. We consider a = 3 , d = (10, 0) and R = 0.5. Numerical 



regions are derived by numerical integration of (81 1, while the numerical roots and the closed 



form plots of the concavity region come from p3\ and (94). When considering the optimality 



of Pr = 1, the numerical roots are obtained by solving (95) and the closed form expression 



through (96). We observe that the conditions derived in Theorems 4.7 and 4.8 predict very well 
the actual concavity regions and regions of optimality of pr = I. In addition, the simple closed 
form expressions are very close to the numerical roots. 

Finally in Fig. [16] we study the maximum rate achievable for the mixed scheme with Pr = I 



for a fixed OP constraint, relative to the DT rate obtained for the same OP. This is equivalent 
to comparing the scenario in which all clusters use DF (pr = 1) with the scenario of using DT 
(Pr = 0). When the ratio of the maximum rates falls below one, then p^ = is better than 
Pr = 1. We consider A^ = lO"'*, d = (10, 0) and a = 4. The rates come from using (55) in (91 ). 



We see that DF provides increasing gains as the relays become more concentrated around their 



sources. However, these gains are smaller than the ones observed in section III-B where only 
the source at the origin had a relay. This shows that the aggregate interference introduced in the 
network by the cooperation can reduce its benefits significantly. 



Finally in Fig. 17 we see the relative reduction of k{R) with respect to kd-p{R) as function 



of (Tin/D for several cases of interesjj In particular we consider the case of Kmix(-R) which 
corresponds to the mixed randomized protocol analyzed above, the case of kdf(R) which 
corresponds to the case in which only the source at the origin has a relay associated (section 



III-B) and kcs{R) (calculated for the case in which all the nodes in the network can have 
an associated relay) which corresponds to the asymptotic value of the lower on the asymptotic 



error probability in ( 1 1 ). Several comments are in order. The value of kdf{R) presents important 
reduction with respect to the case of K]:,t{R) . In fact, for ain = presents a reduction of 100 

^We are making an abuse of notation here, denoting as k{R), the quantities plotted in Fig.llTl In fact they are only the upper 
bounds the corresponding upper and lower bounds derived in the previous sections using the approximations obtained for the 
OP. 
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% with respect to kdt(-R)- In this situation, the relay associated with the source at the origin 
physically collocated with it. This is equivalent to a MISO system with two transmitting antennas. 



The plot and the definition of k{R) in equation (12) suggests that the gain observed is similar in 



nature to a diversity gain of 2 [25 1. As crj„ is increased, we see that the gain diminishes, as the 



second antenna (the relay), cannot decode for free the source message, in order to achieve the 



second term in (14). The behaviour of fi;mix(-R) is similar. However, as in this case all sources in 
the network can have an associated relay, the gains are smaller. This is because, a given source- 
relay-destination triplet has to deal with more interference. The difference between the curves 
corresponding to kdf(-R) and fi;mix(-R) can be interpreted as the penalty of allowing cooperation 
network-wide. We see that, when crj„/-D ^ 0.27, there is no reduction of fi:mix(-R) with respect 
to the baseline of kdtI-R) • This is because in that case the optimal value of Pr is zero, and 
cooperation in not longer beneficial for the network. The value of KcsiR), as obtained from a 
lower bound to the asymptotic error probability, presents a better behaviour with respect to the 
baseline. It is important to notice that when ain is small, the value of K^i^(R) is near optimal. 
However, as a^ increases there is a rapid degradation of K^aiR) with respect to the optimal 
value. 

V. Exploiting CSI at the Relays 

In the above analysis we assumed that each relay decides to activate itself with probability 
Pr independent of all network parameters, which gives rise to a randomized mixed scheme. 
However, a better strategy would be to let the relays decide when it would be convenient to help 
their corresponding sources based on some local knowledge they could have. Given that each 
relay has full CSI of the corresponding link with its source, the relays could determine whether 
they would be able to decode the message from the source or not. In the latter case, the relay can 
communicate his impairment to the destination (at a negligible communication cost) and avoid 
the retransmission of an erroneous message and the corresponding aggregate interference to the 
network. In this way the destination will try to decode the message using only the transmission 
from the source node, which reduces the scheme to DT Clearly, this scheme can not be worse 
than DT. The corresponding OP of this scheme, which uses the full CSI information at the relay. 
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Fig. 14. Outage probability 



,ix(-R, Pr) as a function of Fig. 15. Maximum values of ain/D for concavity of the OP 



Pr for values of ain showing optimality of pr = or pr — 1. 
d = (10, 0), Xs = 10""', 7? = 0.5. Montecarlo simulations 
are obtained by averaging 8 x 10® realizations of the PPP 
using J5I and l|6). Approximations come from using l |55| ) and 
l|56ll in 1I9TI1. 



and for p,. = 1 to be optimal, a — 3, d = (10, 0), R = 0.5. 
Numerical integrals come from numerical integration of (ISTl. 
The numerical roots come and the closed form plots of the 
concavity region come from |93} and j94| l. The numerical 
roots come and the closed form plots of the optimality of 
Pr = 1 come from l|95|l and j96[>. 



can be written a^ 

Pout, Fcsi (R) = inf Pe {BiR, p) n A'^iR, p)} + P© {Adt (R) n A (R, p)} (97) 

pec,|p|<i 

It is important to write the values of the interference power at the relay and destination associated 
with the source at the origin: 





1 U 12 
\Xi - r| " 


i-XiG^s 


h |2 



1{AI{R,P)} 



- + t{AliR,p)} 



\lll..r 



2^{K^rhl^p} 



\xi + ki — r||" \\xi — r\\ 2 \\xi + h — r||^ 

\hkj^ 2'Si {h^^dhlap} 

\xi + ki- d\\" \\xi - d\\^\\xi + ki- d\\' 



, (98) 



(99) 



where Aki (R, p) denotes the event of decoding failure at the relay associated with node Xi'. 



^kAR,p) 



I '^l^X I 

T, \\h-\\(^ 

"■ill '■II 



<T\. 



(100) 



'it should be clear that in this case O does not include 5q and $s does not include the marks Si, because the activation of 
the relay is done with another mechanism. 
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Fig. 16. Maximum rate for the mixed scheme with pr — 1, Fig. 17. Relative reduction of K,nix(i?), i^cs{R) and kdf(-R) 

relative to the same rate of DT for a fixed OP. The mixed (only source at the origin has a relay) with respect to kdt(-R) 

scheme rates are obtained by using i55i in |9l| . As = 10^*. as function of atn/D . R = 0.5, As = 10"'*. d — (10,0), 

d= (10,0), a = 4. a = 4. 



From the previous equation it is observed that the interference power at a given position in space 
depends explicitly on the interference at the points of space where a relay can be positioned. That 
is, in order to obtain the interference power at a given position, we need to know the interference 
power at all relay positions. This could be very difficult to solve, and from a practical point view 
a given relay decision would be linked to all the others relays ones. If a given relay decides not 
to transmit, this would affect the decision of other relays in the network. This means that the use 
of full CSI at each relay leads to a transmission scheme that in essence requires some kind of 
network-wide consensus, which could be impossible to achieve in practice. For this reason we 
will propose a scheme which use partial CSI at the relays. In precise terms, the relay associated 
with the source at the origin, will consider the following set: 



^(0 



\h P 



>e 



(101) 



where ^ is an appropriately chosen threshold. When V (^) happens the relay will try to decode 
the message, because the event V {^) indicates that the fading and path loss gains of its link 
with the source node are strong enough. The outage probability of this scheme can be written 
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as: 



Pout.pcsi(i?)= inf Fe{iB{R,p)UA{R,p))nV{0} + Fe{ADT{R)nV^O} (102) 

pec,|p|<i 

The expressions for the interference powers at the relay and destination associated with the 



source at the origin are given by (98) and (99) with 1 {Ak^ {R, p)} replaced by 1 {T>k^ (0)' 



where V^- (0 is the analogous to (101 ) with the relay associated with source x^. Notice that as 
the events 2)^. (^) are not taking into account the possible impairment of the interference, the 
use of "Dfe, (0 does not guarantee that the relays will decode their sources messages correctly, 
and this scheme does not necessarily degrade to DT as the one that uses full CSI. However, 
although the decisions of the relays based on Pfc. (^) are not fully randomized due to the fact that 
some information is used, they are decoupled. Besides that, they are statistically independent 
because of the independence between marks ki and hs^kr This means that the interferences 
observed at the relay and destination associated for example with the source at the origin are 
statistically equivalent to the randomized mixed scheme of the previous sections with pr chosen 
equal to Pe {V (^)} = e"^"''"". However, the outage performance does not need to be the same 
as the randomized mixed scheme one. This is because, the sets "Dfc. (^) have some information 
about the decoding capabilities of the relays. In fact, it can be shown that this scheme is better 
than the randomized mixed scheme when the interference levels are the same to the relays and 
destinations in both protocols, or what is the same, when pr = Pe {'Dk^ (^)}. In order to prove 
this result notice that because of the independence among the fading coefficients we have: 

Pe {Adt (R) n V' {0} = Pe {Adt (i?)} Pe {V' (0} • (103) 

We can also write: 

Pe {{B{R, p) U A{R, p)) n P (0} = Pe {{B{R, p)nV (0)} 

+ Fe{A{R,p)nV{OnB''{R,p)}. (104) 
It is straightforward to show that: 

Pe {B (R) n P (0} = Pe {B (i?)} Pe {V (0} ■ (105) 

Let us analyze the term Pe {A{R, p) riV (^) fl B^{R, p)}. It can be easily shown that this term 
contains the effect of the partial CSI at the relay, and how the use of that information could 
improve the outage probability performance. We have the following lemma: 
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Lemma 5.1: For every ^ > the following holds: 

¥e{A{R, p)nV (0 n B%R, p)} < Fe{A{R, p)nB%R, p)} Pe {V (0} 



(106) 



Proof: Notice that we can write: 

,{A{R,p)nviOnB'{R,p)} = E. 



P, 



eil-s 



A{R,p)nV{OnB'iR,p)\^,,r 



. (107) 



Conditioned on $s it is clear that A{R, p) fl V (^) is independent of B^{R, p). Using this fact 

we can write: 



, (108) 



jointly with the results from subsection III-A 

Fe{AiR,p)nViOr}B^iR,p)}=E^^^^^[F^^^jAiR,p)nViO\^s,r\Fy{T4[ 

where we have denoted with /| the interference at the destination of the source in the origin 
when the threshold is ^. It can also be shown that: 

T/^ikir f 



i-iPi^ ^ < — 







otherwise 



In this way we can write: 



Fe{AiR,p)nViOnB%R,p)} = e 



-cikir 



xE 



^s,r 






Using the fact that: 

PeM(i?,p)nS^(i?,p)}=E|,^ 
and that for all ^, /| > 0: 



l^AWW^ 



1-e ^^^ \Fv\Tl\i 



iU< 



r/5 



l-e 



iww 



TJ^ II 

q: __ -*- -* T' 1 1 



-l|rho|2 < ]_ — e l-|Pl^ 



1- |p|2^ 

the result in ( 106[ ) follows straightforwardly. 

We have then the following corollary, which follows directly from the above lemma: 
Corollary 5.1: Assume that Pr = Pe {V (^)} = e"^"'''!". Then we have: 

Pout, PCSI (-R) < Pout,mix {R) ■ 



(109) 



(110) 



(111) 



(112) 



(113) 
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This result shows that the scheme using partial CSI has a better performance than the mixed 
randomized scheme analyzed in the previous sections, when the ratio of active relays is the same 
in both (measured by p.,. and Pe {V (^)} respectively) or what is the same, when both schemes 
are subjected to the same interference levels. We should mention, however, that our numerical 
results (not shown) suggest that the performance gains by using partial CSI are not significant 
with respect to the mixed randomized scheme. Moreover, it seems that the performance of the 
scheme using partial CSI is extremely sensitive to the parameter ^, whose optimal value cannot 
be computed in closed form and depends on all the network parameters. 

VI. Summary and Final Remarks 

In this work we presented an analysis of the overall balance of cooperation and interference 
generation in large interference limited wireless networks. The emphasis was put on networks 
where the relays are located in the proximity of their sources. This naturally leads to the use 
of the DF scheme in order to obtain the best possible benefits from cooperation. In fact, the 
principal scheme analyzed in this paper is a mixed randomized one in which DT or DF can be 
used. The choice between DT and DF is done by the corresponding relay associated with each 
source via a randomized decision with probability p,. and without taking into account any piece of 
information the relays could have. This simple procedure, which is mathematically tractable, can 
be thought as a MAC layer at the relays (in a similar fashion as the popular ALOHA protocol), 
with the objective of limiting the interference generation in the network. With this simple model, 
a balance between cooperation and interference generation can be established in the network. 
Surprisingly, for several cases of interest, and for typical operating conditions, the optimal values 
of Pr are or 1, revealing a binary behavior: the best is that all nodes in the network enjoy the 
benefit of cooperation or none at all. There is not an optimal ratio of active aiding relays in the 
network. This means that cooperation could be beneficial or detrimental to all the users in the 
network. Even when cooperation is beneficial to all, the performance improvements are not as 
large in the typical fading relay channel with Gaussian noise. The reason of this comes from the 
fact that, in addition to fading, there is an additional random component given by the network 
nodes positions. Moreover, as discussed in section |V] the use of partial CSI at the relays as 
guidance for their activation decision does not seem to offer additional improvements. And the 
use of full CSI, although surely more beneficial, could be very difficult to implement because 
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some kind of "network consensus" is required. 

Notice that, despite the above results, it should not be implied that cooperation will not 
provide improvements in interference limited scenarios. The fact that relays and destinations 
treat all the interference from the other users as noise severely limits the benefits of cooperation 
for a given source-destination pair. A potential improvement could be obtained using more 
sophisticated cooperative transmission schemes which could take into account the impairments 
generated by the nearby interferers [5] (which introduce by far the most harmful interference). 
In such situation, besides the intrinsic benefits of cooperation, the smart use of the aggregate 
interference introduced in part by the cooperating nodes, could ameliorate its harmful effect on 
the overall network. Also, the study of other cooperative schemes as AF and CF deserves full 
consideration. All these issues, as well as the effect of using several potential relays instead of 
only one, constitute important and interesting future work directions. 

Appendix A 
Interference RVs and their Laplace transforms 

This appendix is presented to simplify the use of the Laplace transform of interference RVs 



throughout this work. For a more detailed analysis see [13| and [14|. Let $ = {{xi,m)} be 
an independently marked HPPP with $ = {xj} an HP?? in M^ and m a vector of marks on a 
subset of M', / > 1. Define the interference RVs: 

Id= ^ fi{d,Xi,mi) /^ = ^ f2{r,Xi,mi) (114) 

where /i and /2 are real valued non negative functions. Then the joint Laplace transform of the 
interference RVs at point ((jJi,u;2) is known to be |13|: 

Cj„Uui,U2) = exp {-A.y^ (1 - E„ [e-'-^M'i^-^"^)-'-^Mr,.,m)^^^ ^A (^^g-. 

Taking ui = or U2 = the single Laplace transforms are obtained. 

Lemma A.l: When fi{d,x,'m) = \hi\'^l(x,d) and f2{r,x,m) = |/i2p/(x,r) with |/iip and 
|/i2p independent exponential RVs with unit mean and /(-, ■) is the simplified path loss function. 



the joint Laplace transform (115) can be written as: 

Cj„iM,uJ2) = e-^4^(-?'''+-^''")+^(---^)], (116) 
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where f{ui,U2) is given by (34), C is given by (35) and T{z) = J^ t^ ^e ^dt is the usual 
Gamma function. 



Proof: In this case m = (|/iip, 1/^2^) so that: 



Ern[e 



-i^lfl{zi,x)-U}2f2{z2,x) 



1 



{loiI{x, zi) + 1){lo21{x, Z2) 



(117) 



Replace (117) in (115) and factorize the integrand as: 

1 1 1 



+ 



{uJil{x,Zi) + l){uj2l{x,Z2) + l) l+M{x,Zi)-^ 1 + M{X,Z2)-^ 



OJl 



102 



(l+J_/(^,^^)-l)(l + J_/(^,^2)-l) 



For the simplified path loss function we use that: 

1 



-dx = Cloi 



2/a 



(118) 



(119) 



l + {ujil{x,r)y 

a result which is known from the direct transmission case [|3|. ■ 

Lemma A.2: Suppose the marks of the HPPP are m = (|^ip, \h2\'^,e,k), with |/iip and |/i2p 
unit mean independent exponential RVs, e a Bernoulli RV with success probability pr, and k 
a RV on M^. Let fi{d,x,m) = {\hi\'^ + e\h2\'^) l{x + Tk,d) with /(■,■) the simplified path loss 
function and r G [0, 1]. Then the Laplace transform is: 



£/^(wi) = exp <^ -XsCu/ 1 + 



2/a 



2pr 

a 



(120) 



Proof: Take u;2 = in (115) and compute the expectation with respect to the marks: 



^m[e 



-uiifi{d,x,rn)~\ 



Pr 



[1 + coil{x + rk, 



jdFk + 



1 ~Pr 



1 + ujil{x + rk, d) 



dFk. 



(121) 



Replacing ( |121[ ) in ( |115[ ) and interchanging the integration order we find ( |115[ ) yields: 

1 



£/^(wi) =exp<^ -XsPr / ( 



1 



[l + ujil{x + Tk,d)f 

Xs{l -Pr) 



dx I dFk 
1 



1 + {uil{x + Tk,d)y^ 



dx dFk 



(122) 
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When the integrals with respect to x are computed the result does not depend on k so the 
distribution of k does not affect the final result. For the first integral we have: 

1 ^dx = 27ra;2/° / —— — -tdt (123) 

^ + 2./ (TTT^*). (124) 

= u,'i'° (l + ^) C. (125) 



For the last step we integrate by parts and C is defined on (35 1. The second integral is (119). 



Appendix B 
Proofs of section [m] 



A. Proof of lemma 3.1 



We need to prove that p = minimizes the following function: 



|2 



1-^ 1- 



J[P) = T\ 1 |2^2M + ' (^2^) 

(1 _ |p|2)2/a /i2-/ii 



where /^i and yU2 are given by ( [22] ) and ( |23| ). Since the first term is minimized when p = we 
focus on the second one. Define: 

a = lsd + lrd, b = [{Isd - Irdf + ^lsJrd\p\'^] ■ (127) 

SO that we can write /ii = | (a — 6) and ^2 = | (a + 6). Notice that < 6 < a, that p appears 
only in b and that b is increasing in \p\. Therefore it suffices to show that the auxiliary function: 

is increasing in b to complete the proof. We differentiate to obtain: 

dg.,. 1 ia-b)-^a-H)-ia + b)-^a + ^) 

db^'^~2'--. b^ ■ ^ ^ 

In order to conclude the proof we need to check that dg/db{b) > 0. That is: 

t{b) > u{b), (130) 

where 

t{b) = {a + b)^ (a- -bj , u{b) = {a- b)^ (a + -bj . (131) 

It can be readily proved that t(0) = u{0). Through standard calculations, it can also be proved 
that when a > 2 and b < a, dt/db{b) > du/db{b). These results, permit us to conclude that 



t{b) > u{b), or what is the same, the second term in (126) is minimized when p = 
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B. Proof of lemma \3.3 



For any u > 0, I) > and a > A define the auxiliary function: 



{u):- 



gD2(l-«2) _ ^c 



(132) 



Now we show that by replacing 



9[u) 



.DHl-u^) 



E 

n=0 



n\ 



(133) 



in (132) by griu) := 1 + ^'^(l - u^) + ^(1 - u^f we obtain a lower bound of (132). When 



u G [0, 1] the series in (133) is monotonously increasing so in fact griu) < g{u). For u > 1 we 
have: 

rlr, _ 

-2D\g{u) > -2D\. (134) 



dg 



du 



[u] 



Now integrate in the interval [1, u] on both sides of ( 134 ), add g{l) = 1 and apply the fundamental 
theorem of calculus to show that g{u) > 1 + _D^(1 — v?). Repeating the procedure of multiplying 
by —2D'^u on both sides and integrating we find that if m > 1 then g(u) < gxiu)- Using this 
two inequalities we find that for all li > we have: 

griu) - m" 



u)> 



l-u° 



(135) 



0(m) and its first derivative have an avoidable singularity at m = 1, so using L'Hopital's rule we 
compute the first order Taylor expansion of (f)(u) at m = 1, 0r(w) = —Mu + N, where: 



N = 1 + D^ + 



2D^ 



M = D' 



a 



a 



-il-D' 



are both positive for a > 2 and -D > 0. Now we show that m > 0: 

gT{u) - u°' 



1 -M° 



> -Mu + A^. 



(136) 



(137) 



This is equivalent to: 



gT{u) - u 



-Mu + N){1 - u'') 



> < u< 1 



< M > 1, 



(138) 



which means that we need to study the sign of: 

p(m, a) = -Mu''^^ + (N- 1)m" + — M^ - (D^ + DV + Mu + D , ^ 

2 V 2 a 



(139) 
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First we show that p(u, a) is increasing in a when < m < 1 and decreasing in a when m > 1; 
to do this show that ^ > when < m < 1 and ^ < when u > 1. Differentiate p with 

da — — da ^ 

respect to a to find that: 



dp_D^ 
da a 



2 



hog{u)u"{l-u)a'^ + 2 {l-u)D^ + u [(log(u") - 1) u" + 1]| . (140) 



Since we are interested in the sign of the derivative, we can neglect the factor D"^ jo? . Notice 
that: 



\-u)D^^u{ (141) 



>u ifO<M<l 

<u if M > 1, 

and in addition, [(log(M") — 1) m" + 1] > when n > 0, attaining a minimum at m = 1. This 
can be proven by differentiating: 

iKM!^!^flli^l±i!=a."-MogK), (142) 

an 

which implies that this function decreases when < m < 1 and increases when m > 1. Using 

this fact, ( 141| ) and defining: 

/(u, a) = log(M)u"(l - u)a^ + 2u [(logK) - 1) m° + 1] (143) 

we can bound: 

a^rfpj>/Ka) ifO<M<l 

^'^" [</(«,«) if«>l. 

When < M < 1 we have that /(m, a) is increasing in a and when m > 1 it is decreasing. To 
see this, differentiate: 

-^ = an" log(u) [2(1 - u) + log(u)(a(l - u) + 2u)\ (145) 

da 

an observe that for m > 0: 

2(1 -u) + log(M)(a(l -u) + 2u) < 2(1 - M + log(M)) < (146) 

so the sign of df /da is the opposite of log(M) for each pair (u, a). This means that: 

> f(u,2) ifO<M<l 
fiu,a){- - (147) 

< f{u,2) ifM>l. 
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Finally, continue differentiating to show that /(m,2) > if < m < 1 and f{u,2) < when 
u > 1. This means we proved that: 



a^dp J >/(«,«)>/(«, 2) >0 ifO<M<l 



(148) 



< f{u, a) < f{u, 2) < if M > 1, 
which concludes the on the monotonicity properties of p{u, a) with respect to a. This implies 



that it suffices to show that condition ( 138 1 is met when a = 4. Notice that in that case we have 



M = f (1 + D), N = 1 + D + D"^, which imply that p{u) can be factored as: 



p{u) 



-^{D + D'^)uiu-lf{u + l). 



(149) 



We conclude that (137) is true by observing that (149) satisfies condition (138). To conclude 



the proof of the lemma use (135), (137), take 



r — 



u 



D 



D = {XsSy/'D, 
and multiply both sides of the resulting inequality by e"'^^'^^ . 



(150) 
(151) 



C. Complement to the proof of theorem \3.3 
To find ( [55] ) start by writing: 

EJIIr- 



l _ikir 

"r — d\\e ^'^«" dr. 



27ra_ 
Now take x = r — d, change to polar coordinates to obtain: 

E,.|||r-<i|||=y^ (-) e-^/„ (^ j *a 

= CrQ2.o { — ,0 



(152) 

(153) 
(154) 



In the first step we used the definition of the modified Bessel function. For the actual value of 



Q2,o(m,0) use (91) from [31 1 which states that: 

" POO 

Q2fl{b,ac) = a / Qi{b,ax)dx + cQi{b,ac) 

-J c 

Taking c = 0: 



Q2,o(^)0) = a Qi{b,ax)dx. 
Jo 



(155) 



(156) 
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and using (60) of pTJ which states that: 

) 

Qi{b, ax)dx ■■ 



27r _b^ 
e 4 



4a 



(^^ + 2)/o(^)+6^/i(^ 



(157) 



we obtain the desired result. 



Finally for the upper bound on E,. [\\r — d\\\ start from (152), change to polar coordinates and 



apply lemma D.l To solve the resulting integral use integration by parts and that: 

2sin( - J -1 J d^ = 2(4-7r) 

uf{u)du = — :^e~2^ 



27r 



u 



u f{u)du = —^r-e 2,72 -)- 



a 



2-K 



I erf 



Stt \V^a 



u 



(158) 
(159) 
(160) 



(160) is obtained integrating by parts. 



Appendix C 
Proofs of section [IV] 



A. Proof of theorem 4.2 



Rewrite Pout(-R,Pr) in y\\ as: 

Pout(i?,Pr-) ^ 



l-exp<^-A,5D^ 1 + 



2pr 

a 



+ 



PrD' 



D'^-Wr -d\ 



exp <^ ~\,5D^ 1 + 



2pr 



a 



exp<^-A,5||r-rff 1 + 



2pr 



a 



. (161) 



It is immediate to show that the first term in (161) is concave in pr. Hence, only the concavity 
of the second term is to be proved: 



h{Pr 



PrD' 



D" - \\r -d\\' 



exp <^ -A,(5D2 1 + 



2pr 



a 



-exp{ -XJ\\r -df 1 + 



2pr 



a 



(162) 



The second derivative dh'^ / dpl{pr) of h{pr) can be written as: 



d'^h 
dpi 



4\JD' 



a{D°' - \\r-d\\°') 



Pr 



a 



\s5\\r - df) \\r - rff exp <^ -\s5\\r - df ( 1 + 



2pr 



a 



(l - ^KSD^) D^exp l-XJD^ (^ + ^) | 



(163) 
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Let us analyze the value of the function j^{pr) when pr = 0. We obtain: 

W' = a(D'-\\r- dr ) [II'- - ''I' '^-P {-^-^ll'- -m-D'e.p {->.,W'}] , (164) 
Consider the case u = ^^''^2" < 1. It is easy to verify that (Ph/du'^{0) < if: 

111 11 

XsSD^ < , ue\0,l]. (165) 

1 — u 

It can be shown that —j^ is decreasing in [0, 1] with 

\nu Inu 
lim = +00, lim = 1. (166) 

■u-s-O 1 — U "^>1 1 — U 

This means that when u G [0, 1], ^(0) < if XgSD"^ < 1. The case u > 1 can be handled in 
the same way to obtain that when m > 1, ^(0) < if \s5\\r — rfp < 1. Finally, the case u = 1 



can be treated with appropriate continuity arguments. Summarizing, we have that 4-|(0) < if 



dpi 

dpi 

A,5max{DM|r-c/f } < 1. (167) 



As h{pr) is continuous function, the fact that ^(0) < permit us to assure that there exists an 
interval [0, (3] with /3 > where 4-|(Pr) < if p.,. G [0, /?]. In order to guarantee the concavity of 



h{pr) for Pr G [0, 1] we need to find conditions for /3 > 1. We need to consider (163). Consider 

d^hi 



again m < 1. In order to find the points p*. where |-|(p*) = we need to solve: 

gipl) = exp l-Xs6D' (l + ^ Vl - u)\, (168) 



a 
where 

It can be easily shown that g{pr) is an increasing function of Pr is n < 1. Note that g{pr) has 
a vertical asymptote at j-frp ^^'^ ^^^° ^^ x s\\r-dp - ^^ ^ ' ^ " ^ is decreasing in pr 



when u < 1, it is not difficult to see that there exist only two values p* such that (168) is 
satisfied, and such points satisfy: 

(1) " p(2) > ^ (170) 



(2) 
From the condition (167) and the fact that a > 2 it is clear that pr is of not use for us. We 



just need to concentrate on proving that pi > 1- As g{pr) is an increasing function of Pr, we 
just need to prove: 

^(1) < exp !^-\s6D^ (1 + ^) (1 - u)\. (171) 
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Consider the existence of a value K > 1 such that: 

. XsS\\r~df 



1- 



XJD2 



< K. 



It is straightforward to obtain: 



\JD^ < 



a{K-l) 
K -u 



From (171 1 and (172) we obtain also the following condition: 



K5D^ < -7^-^ 



1 



(1 + !)(!-«)' ^^^' 



(172) 



(173) 



(174) 



Combining (173) and (174) we obtain the following condition which guarantee us that (171 ) be 
satisfied: 

K^D"^ < (3% (175) 



where /3* is defined as /3* = infu<i/3(n) and 

'a{K-l] 



/3(n) 



max mm 

1<K<^ 



In (Ku) 



(176) 



K-u {l + i)il-u) 
Notice that (3* depends only on a and can be optimized numerically. This numerical optimization 
suggests that the value of /3* = 0.38 for « = 2 is a lower bound for all the values of (3* for 
other values of a > 2. In this way the can conclude that h{pr) is concave as a function of pr if 
XgSD"^ < 0.38. The case n > 1 is treated in the same way. The conclusion of the lemma follows 
straightforwardly. 



B. Proof of lemma 4.2 



For the first item consider that: 






a 



In order to be negative we need: 



|rf-D"^'-|"^-^l"' 



D" - \\r -(i||" 



Irlr < £>"- 



L)2_ ||r-d| 



(177) 



(178) 



Defining x 



D" - ||r-rf||" 

^'''^2 ^ we obtain that in order to ^jf^{Pr) to be negative the following two 
conditions have to be satisfied in r G M^: 

1 — X 



\r - d\\ = xD , ||r|| < D '-^, for some s > 0. 

1 — X2 



(179) 
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This conclude the proof of the first item. For the second we have: 



m(0) = XsSD^' 



m(l) = A,5 1 + 



a 



D' 



\r\'-D' 



D^-\\r-d\\' 



(180) 
(181) 



D" - \\r - d\\' 

In order to see when p^ = 1 is optimal we need to consider m(0) — m(pr) > 0. It is easy to see 
that this happens when: 



D' 



I — X 



1 



2D' 
> 

a / a 



1 



1 — X2 \ a / a ,^ ^ 

Rearranging terms and using terms, we have: 



\r - df = xD^, for some x > 0. (182) 



Irlr <D' 



1 — X 



_l-xi a + 2 
This concludes the proof of the lemma. 



, \\r — d\\'^ = xD'^, for some x > 0. 



(183) 



C. Proof of theorem 4.6 



First we have to evaluate the OP conditioned on the position of the NN of the source at the 



origin, that is, evaluate (90) without the expectation with respect to r. Notice that in that case 
the expression of the OP in terms of the Laplace transforms of the interference RVs is the same 
as in the fixed relay case, and therefore ( [8T] ), which we repeat here for comfort, is also valid: 



IPout,mix(^,Pr) = (1 - Pr) [1 " -^/d {T/lsd)] + Pr 



D 



||r-(i||° D" - ||r -(i||" 

+ Pr[l-CiAT/lsr)]. (184) 

When the Laplace transform is evaluated we have to take into account that the nearest neighbor 



of each cluster is also random. However considering the approximated interference RVs (64) 



and (63) and lemma A. 2 in appendix [A| the Laplace transform in this case are also: 



£/^(wi) = exp <^ -XsCuj/ 1 + 



2/a 



2pr 



a 



(185) 



where 6 is given by pm and C is defined as in (35). Now we have to average (184) with 



respect to the position of the origin's nearest neighbor. It is clear that the first term of (36) is 
not dependent on r. In the second term taking u = \\r — d\\/D we can write: 



D° - ||r-d||' 



\r-drCjAT/ls,) 
/}" - ||r-d|h 



-AsAD2 



1 -U" 



(186) 
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Using first that e^-^^'(i-"') = l + XsAD\l-u^)+0{{XsAD\l-u^)y) as {XsAD^{l-u^)) -^ 
and then that for u > 0: 



u 



1 -M° 



> 1 



1 u, 



we find that: 

gA.AD2(l-„2) 



U 



> 1 + XsAD^ 



■u^ 



1 + 



a 



1 ]u 



O (XsAD 



2\2 



'l-U 



2^2 



M^ 



(187) 



(188) 



as XsAD'^{l — u^) — )■ 0. Taking expectation with respect to m = ||r — d\\/D we obtain: 



E. 



gA,AD2(l-M2) _ 



U 



w- 



> 1 + XsAD' 



1 + 



a 



1 Er [n] 



+ 0{{X,AD^f). (189) 



The third term can be evaluated in a straightforward manner as in lemma 3.3 to yield: 

E, [Cl^ (T/lsr)] ^'^ 



TxXin + XgA 



(190) 



D. Proof of theorem \4.7\ 

Define i'{pr) '■= XsA{pr)D'^ to write the OP as: 

a(z/(p,,) - uiS))) 



■rout,mix(-n-) ^ 



1 - 



[1 



-v{pr) 



2z/(0) 

a{u{j>r) - ^^(0)) _^(p,) 

2z/(0) 



] + 



a{v{pr) - ym 



1 + 



V{pr 



1 + z/(pr 



2z/(0) 



2 — «„ n, 



nXinD^ + iy{pr 

. (191) 



It is clear that analyzing the concavity with respect to ^{pr) is the same as analyzing it with 
respect to pr. Our goal is to differentiate twice with respect to u to find conditions for the OP 
to be a concave function of pr. We have found the exact value of E^ [||r — d\\] in terms of the 



Nuttall Q2fi function (55) but this expression is not easy to handle. For this reason we would 



like to use a suitable bound for E,. [\\r — d\\], that is, we have to differentiate twice and then 
bound the expectation accordingly. For this reason we first find conditions that guarantee that an 
upper bound on the expectation will yield an upper bound of the second derivative of the OP. 
Let us define the auxiliary function: 

(j){u{pr),u) := -e-'^P^\u{pr) - z/(0)) (1 + u{pr)u) (192) 



which is similar to the third term of ( 191 1, and study its behaviour with u. Differentiate 0(z/, u) 
twice with respect to u and then with respect to u to find: 

d f d^c 

du 



(H) = e-'' [-u' + (4 + umu - 2(1 + urn] 



(193) 
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where Xs6 = iy{0) <i'{pr) <'^(1) = ^s^ {^ + "^ ■ Notice that the derivative does not depend on 
u. It suffices to study the values of the function at the borders of the interval of u to see how 
the derivative behaves with u. When p^ = we have v = \s5 and we find that the derivative 
is negative as long as XgSD"^ < 1. When pr = I then u = X^S (l + ^) and we find that the 
derivative is negative if: 



z/(0) = XJD^ < 



A + a 



l-(l + 4/a + 8/a 



2\l/2 



(194) 



2(l + 2/a) 

It is easy to show that the expression on the right is increasing in a which means that the 
function attains a minimum at a; = 2. To do this show that the numerator is increasing in a and 
the denominator is decreasing in a, and that they are both positive. Therefore setting a = 2 in 



(194) gives a sufficient condition for the second derivative to be negative in u in all the interval: 

3- v^ 



1/(0) = Xs6D^ < 



0.38. 



(195) 



Therefore if condition ( 195 1 is met it is easy to show that since a > 2 we can differentiate ( 191 1 



twice with respect to z/, then use ( |56| ) and study the conditions for the concavity of the OP with 
respect to p,.. Now do this and reorder the terms to find: 



t^'Pout.mix ^ anX^nP^ uiO)+TlXinD' 
rfz/2 - 

with: 



-iy{Pr) 



Z/(0) (l^iPr) + TrXinD 



2^3 



HO) 



[A2{sMpr)^ + A,{sMpr) + Ao{s)] (196) 



A,(s] 



Ms) - 

Aois) = 2 



1+1 



1- |)7(s)s, 
a 



a 



7(s)s [4 + z/(0)] 



7(s)s [l + z/(0)] + z/(0). 



(197) 
(198) 
(199) 



and s = ain/D. Now apply the bound: 

anXinD^ u{0) + nX^^D^ 



< 



a 



2as^ 



U{0) {U{pr) + 7rA,„D2)3 - zy(0)7rAi„D2 ^,(0) ' 

and use that 2afj,^ = (vrAm)^^ which is straightforward from the fact that i/(pr) > ^^(0) 
Xs6D^ > 0. Then factor out the term e'^^P'^yiyiO) to obtain: 



(200) 



" ■"■ out,mix ^ c- 



(iz/2 (p^ 



KO) 



[2ae'^(^'-)s2 - Ao{s) - A,{s)u{pr) - A2{s)u\pr)] 



(201) 
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Now using that a > 2 we upper bound: 



-Ai 



1 



a 



7S 



[4 + z/(0)] < 4 + z/(0). 



(202) 



Now using (202) in (201 ) we find that we can find sufficient conditions on s for the concavity 
of the OP by studying the sign of the function: 



PrXuiPr), S) := 2«e"(P'-)s2 _ Ao + (4 + u{0))u{pr) - A2u\pr). 



(203) 



To do this, differentiate (203) with respect to z/ to find that pc{iy,s) is increasing in u for all 



s > whenever (195) is met. This means that a sufficient condition for the concavity of the OP 
is that pc ('^(l), s) < 0. Now we can rewrite pc explicitly as a function of s, evaluate at z/(l) 
and remembering that z/(l) = z^(0)(l + 2/ a), upper bound it: 

p,(z/(l), s) = 2ae-^^h^ + (a - 2) ("^ + z/(0) + l\ ^{s)s + Co(z/(l)), (204) 



< 2ae'^(i)s2 + (a - 2)e''^^h{s)s + Co(z/(l)), 



(205) 



where Cq{v{1)) = -//^(l) + (4 + z/(0))z/(l) - (2 + 3z/(0)). First notice that Co(z/(l)) < for 
a > 2 and as long as z/(0) < 1/2. To show this notice that Co{iy{l)) can be written as: 

1 



Co Ml)) 



a' 



[(z/(0) - 2)a^ + 2z/(0)(4 - z/(0))a - 4z/(0)2] 



(206) 



which we want to prove is negative for a > 2. Since z/(0) satisfies (195) we deduce that this 
function has a single positive root in a, which is: 

"4-1/(0) v^z/(0)2 - 4z/(0) + ! 



ao = '^(O) 



(207) 



2-z/(0) 2-z/(0) 

and we want to prove that ao < 2. Differentiate the function between brackets in ( |207| ) to show 
that it is increasing in z/(0), and therefore, Oq is increasing in z/(0). Taking t> = 1/2 shows that 



in that case ao = 2 so we get that Co (^^(1)) < when (195) is met. Since the other two terms 



of (205) are positive for all z/ > it is clear that there will be an upper bound on s for the 
OP to be a concave function of Pr- Finally, a sufficient condition for the OP to be a concave 



function is finding the positive root in s of the upper bound in (205): 



.AsAD^ 



{2as^ + {a- 2)7(s)s) - {Xs^D^f + (4 + XJD^)XsAD^ - (2 + 3Xs6D^) = 0. (208) 
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Notice that since 7(s) is transcendental function of s we cannot find the roots explicitly. However 
we can upper bound the value of 7(5) in the following way. For each a it is clear that the 



maximum value of s is attained when \s6D = 0, and in that case (208) becomes: 



2as^ + (a-2)7(s)s-2 = 0. 



(209) 



It is easy to show that the roots in terms of s are decreasing with a, and therefore for a = 2 the 
root is the largest over all a > 2 and XgdV^. In that case we have that (209) is independent of 



7(5) and the value of s is s = 2 ^/^. Therefore, solving numerically we find that for all a > 2 



and As(5-D^ > we can take 7(5) < 1/2. In that case (208) is a second degree polynomial in s 



and its roots can be found in closed form. To find (94) take the positive root of this polynomial 
and find a lower bound by discarding the square of the linear coefficient inside the square root 
of the expression. 



E. Proof of theorem 4.8 



We want to find conditions on ai„ that guarantee that Pout,DF(Pr = 1) < IPout,DF(Pr = 0). 
Setting Pr = and p^ = 1 we can write: 



■It out.mix (Pr — Uj — i 6 



-XsSD^ 



■It out,mix (,Pr — 1) — 1 



AsA 



A^A + vrAj 



A, AD' 1 



a. 



aD 



-ErWlr 



-AsAD2 



(210) 
(211) 



An explicit for the values for which Pout,mix(Pr = 1) < IPout,mix(Pr = 0) cannot be found. Therefore 
we find a sufficient condition by bounding ( |210| ) and (211). First we use (56) and: 



A.A 



AsA + nXin 7r\n 



A.A , , 9 
< ^ = 2A.,Aai 



(212) 



in (211 ) to upper bound: 

Pout,mix(Pr = !)<! + 2XsAal 



1 + X.AD' 



1 + (1 + 7 sHt 

a \ D 



^-x.AD\ (213) 



Therefore a sufficient condition for pr = 1 io be optimal is obtained by asking that the right 
side of ( |213| ) to be smaller than (210). Writing this condition an reordering yields that Pr = I 
will be optimal when: 



2XsAD\s' + XsADU 1 



a 



-XsAD- 



7(s)s + e' 



-XsSD'^ 



-XsAD^ 



1 + -XsAD' < 0, (214) 
a 
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with s = ain/D. Now upper bound the terms that do not depend on s as: 

'-XJD^ - 1 



g-A,5D2 _ ^-XsAD'^ 



1 + -XsAD^ 1 < \ e 



a 



su±^ -AsAD^ ' 2 



a^ 



(215) 



using that when condition (195) is met then e" < l + u+ ^u . This turns the sufficient condition 



(214) into: 



2a {a + 2) e^^^^' s^ + {a^ - 4) -f{s)s + 4 ('^KSD^ - 1 ) < 0, (216) 



Again by setting equality in (216) we can find the maximum value of s so that Pr = 1 i^ optimal. 



However this equation does not have an explicit solution and must be solved numerically. To 
work around this problem and obtain a closed form expression, observe that for each a as XgAD^ 



grows the maximum value of s that satisfies (216) decreases so that the maximum value of s is 
attained when Xs5D'^ = 0, which is also intuitively right. Therefore setting XgdV^ = in (216) 



yields an upper bound for the condition for each a: 

2a {a + 2) s^ + (a^ - 4) 7(5)5 - 4 < 0. 



(217) 



It is easy to see that as a grows the maximum value of s that satisfies the equality in (217) 



decreases and therefore setting a = 2 yields an overall upper bound for s for all a > 2 and 
Xs5D^ > 0. Doing this we obtain that s < 1/2 which yields 7(5) < 2/5. Upper bounding 7 in 



(216) turns the condition the equation into a second degree polynomial in s that can be solved in 



closed form. Once more to find (96) take the positive root of this polynomial and lower bound 
it by neglecting the square of the linear coefficient inside the square root of the expression. 



Appendix D 

Auxiliary results 

A. CCDF of a quadratic form of complex circularly symmetric Gaussian RVs 

We want to calculate the distribution of: V := \hsd\'^lsd + \hrd\'^lrd + 2y/lsdlrd^iphsdKd)^ which 
can be factorized as: 



V^ = g^ 



hd Vhd^rdP 

yhdlrdP* Ird 



q ■= q^Qq, 



(218) 



where q 



h h 



is a zero-mean complex circularly symmetric Gaussian vector with 



identity covariance matrix. Since Q is positive definite we can diagonalize it to obtain V 
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q^ PDP^q, where P is a unitary matrix. Defining w 



Wi W2 



Pq we have that 



V 



w^Dw = Dii 



Wi\ 



-D2,2|w^2p, where Dii and D2,2 are the elements on the main 



diagonal of D. Since P is unitary, w = P^q has the same distribution as q so {|tfip, |u^2p} are 
i.i.d. exponential RVs with unit mean. This means that Di^i\wi\'^ and -D2,2|^^2p are independent 
exponential RVs with means Dii and D2,2 respectively, where: 



^1,1 



D 



2,2 



1 /9 

(^sd + Ird) + ((^sd — ^rd) + 4:lsJrd\p\ ) 
(hd + Ud) — {{hd — Ird) + ^hdlrd\p\ ) 



(219) 
(220) 



are the eigenvalues of Q which are on the main diagonal of Q. In Appendix D-B taking V 

X + Y , jjii = Di^i and /i2 = -D2,2 we obtain: 



Fviz) 



M2-M1 



/il 7^ /U2 

1 - (1 + 2;//ii)e-^/'^i /ii = /i2. 



(221) 



Observe that /^i = 1x2 if and only if Isd = U-d and p = 0. 



5. CDF 0/ f/ze sum of independent exponential RVs 

Let {X, Y} be independent exponential RVs with means pi and ^2 respectively. The distri 
bution of their sum is: 

. If 111 = H2- 

Fx+y{z) = 

Jo 



If Ail 7^ P2- 



Finally: 



' fx{.x) { / fY{y)dy\ dx 



Fx+y(z) 



fx{x) ( / fY{y)dy ) dx 



1- 



^2e-''^^ - /iie-^/'^i 



'X+Y\ 



P2 - Pi 

M2-M1 



Pi 7^ P2 



1 - (1 + z/fii)e-'/''^ pi = p2. 



(222) 
(223) 



(224) 
(225) 

(226) 
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C. An additional bound 

Lemma D.l: Let r and d be fixed vector in M^ and let 9 E [0, 2n) be tlie angle between r 
and d. Then the following bound is valid: 

\\d-r\\ < \\\d\\- ||r|||+2min(||rf||,||r||)|sin(^/2)|. (227) 
Proof: Write: 

r — d = Ui + Vi (228) 

r-d = U2 + V2, (229) 
where: 

ui = \\r\\ (^ - ^1 ^1 = ^ (Ikll - IMII) (230) 



«2 = IMII ( ^ - A I V2 = Tr^ i\\r\\ - \\d\\) . (231) 



\r\ 
Use the Pythagorean theorem to show that: 

||mi|| =2||r||sin(e/2) (232) 

lluall =2||rf||sin(^/2) (233) 

ll^i|| = 11^211 = |||r||-||rf|||. (234) 



Finally apply the triangle inequality to (228) and (229) to find: 

\\r-d\\ < \\\r\\ - ||rf|||+2||r||sin(^/2) (235) 

\\r-d\\ < \\\r\\ - ||d|||+2||d||sin(^/2). (236) 

To conclude we can take the minimum between ||r|| and ||rf|| from the previous equations. ■ 
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